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ABSTRACT 


Frequently  the  applied  statistician  is  confronted  with  the  problem 
of  investigating  the  distribution  of  experimental  data  sampled  from  a 
continuous  two-dimensional  distribution.  A  prime  example  in  the  area  of 
weapons  system  analysis  is  the  investigation  of  the  distribution  of  hit 
points  of  missiles  or  projectiles.  Because  of  the  desirable  characteristics 
of  the  bivariate  normal  distribution  and  because  of  the  frequency  of 
occurrence  of  these  characteristics  in  physical  phenomena,  the  applied 
statistician  is  often  interested  in  determining  if  the  sample  distribution 
of  the  experimental  data  can  be  approximated  by  a  bivariate  normal  distribution. 
The  Bl-CHI  (Bivariate  Chi-Square  Test)  computer  program  described  in  this 
report  is  for  assessing  the  agreement  of  a  continuous  two-dimensional  sample 
distribution  with  a  parent  bivariate  normal  distribution.  The  program 
computes  the  Chi-square  statistic  for  testing  the  null  hypothesis  "The  random 
sample  is  from  a  bivariate  normal  parent  population."  The  minimum  and 
maximum  sample  values  of  each  random  variable  are  determined,  and  the  sample 
estimates  of  the  five  parameters  (two  means,  two  variances,  and  the 
correlation  coefficient)  of  the  bivariate  normal  distribution  are  computed. 

Also,  the  sample  estimates  of  the  two  regression  lines  (Xg  on  Xj  and  Xj  on  Xg) 
are  computed.  An  optional  feature  allows  additional  output  of  a  plot  of  the 
sample  data  and  the  contour  ellipses. 

The  program  is  coded  in  FORTRAN  IV  for  the  IBM  7030  (STRETCH)  computer. 


iii 


FOREWORD 

The  formulation  of  the  BI-CHI  computer  program  was  performed  in 
the  Mathematical  Statistics  Branch  of  the  Operations  Research  Division. 
The  coding  of  the  program  was  performed  in  the  Operations  Sciences  Branch 
of  the  Computer  Programming  Division.  The  project  was  supported  under 
Foundational  Research  Project  No.  29Y,  "Computer  Programs  for  Statistical 
Analyses ." 

The  original  program,  used  for  exploratory  computations,  was  coded 
by  Mr.  Rick  Ball;  the  final  program  as  documented  in  this  report  was 
coded  by  Mrs.  Jo  Ann  Brown.  Miss  Nancy  Dean  coded  the  computer  routine 
which  generated  the  c* -values  tabulated  within  the  BI-CHI  program. 

The  date  of  completion  was  31  October  1967. 
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I.  INTRODUCTION 


BI-CHI  (Bivariate  Chi-Square  Test)  Is  a  computer  program  for  the 
Chi-Square  Goodness  of  Fit  Test  for  a  bivariate  normal  distribution. 

The  program  performs  a  statistical  analysis  of  the  discrepancy 
between  sampled  data  from  a  continuous  two-dimensional  distribution 
and  the  hypothesized  bivariate  normal  distribution.  This  discrepancy 
is  measured  by  classifying  the  n  randomly  sampled  (Xj t jXg,) -observations 
into  k  mutually  exclusive  and  exhaustive  categories  (areas  between 
coaxial  contour  ellipses)  and  comparing  the  observed  frequency  with 
cne  theoretical  frequency  of  each  of  the  k  areas.  The  formulation  of 
the  Goodness  of  Fit  Test  is  described  in  Bates  (1966)  and  therefore 
is  not  repeated  in  its  entirety  in  this  report. 

The  joint  probability  density  function  (pdf)  of  two  jointly 
normally  distributed  random  variables  xx  and  is 

f(x,  ,xr)  -  K  'e’^i-  «  <  Xj  <  ®,  -  ®  <  Xj,  <  <»',  (1) 


where 


Q  = 


TT? T 


-2p 


‘i’Ux 


*a  ’Us 

+ 

.2 

CT8  , 

a3  1 

(2) 


A  horizontal  plane  parallel  to  the  x^ -plane  intersect)  the 
surface  of  equation  (1)  in  the  shape  of  an  ellipse  (contour  ellipse). 
The  equation  of  the  general  contour  ellipse  is 
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where  c  is  a  constant.  The  contour  ellipses  are  centered  at  (u1  ,(i2)  , 
and  their  general  shape  is  determined  by  and  p.  By  performing 

an  orthogonal  transformation,  the  quadratic  form  (Q)  can  be  3hown  to 
be  distributed  as  y2  with  two  degrees  of  freedom.  That  is, 

Q~X2<2)-  (4) 

Therefore,  the  probability  that  (Xj.Xj)  falls  in  the  contour  ellipse 
of  equation  (3)  is 


PrXC<2)  <  c2]  =  p, 


(5) 


where  c2  is  defined  by  p  =  I  Pf(ya)dy2,  where  f(ya)  is  the  pdf  of 

J  o 

\2(2).  Additionally,  the  probability  that  (x1 .Xg)  falls  inside 
the  p"  contour  ellipse  and  outside  the  p'  contour  ellipse  is 


prCp»  «  X: 


(2) 


c'2..]  =  p"-p' 


0  p1  <  p"  s  1, 


(6) 


The  classification  of  the  (x; , ^ ( ) -observations  is  then 
performed  by  evaluating  Q,  (i=l, 2,  •• ‘,n)  .  If  c’  *;  Q,  <  ca 

”3-1  1  t3 

( j  l  ,2,  •  •  • ,k) ,  the  ith  observation  falls  within  the  jth  contour 
ellipse  and  outside  the  ( j - 1) C ^  contour  ellipse.  Denoting  the 
observed  number  of  (x  falling  within  the  jth  contour  ellipse 

and  outside  the  ( J - 1) C ^  contour  ellipse  by  0,  and  the  expected 
number  by  E  •  n(p  . -n,  ),  the  statistic 


k 

X  r  (O.-E.)VE.  (7) 

j=i 


is  approximate- ly  distributed  as  y?(k-l>  . 
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In  practice,  however,  the  population  parameters  of  the  hypothesized 
bivariate  normal  distribution  are  unknown.  That  is,  the  population 
is  not  completely  specified  under  the  null  hypothesis;  only  the  type 
of  the  parent  population  is  known  or  assumed.  In  such  case,  the 
population  parameters  are  estimated  from  the  sample  data;  Q(  is 
replaced  with  Q, ;  and  equation  (7)  becomes 

k 

X"  =  I  (0, -E ,  (8) 

3  =  1 

A 

where  E,  denotes  the  use  of  sample  statistics  to  estimate  population 
parameters.  The  degrees  of  freedom  of  X3  is  reduced  by  one  for  each 
parameter  estimated.  Therefore,  when  all  five  parameters  of  the 
bivariate  normal  distribution  are  estimated,  X3  is  approximately 
distributed  as  xr(k-6). 

BI-CHI  computes  the  five  estimators  of  the  population  parameters, 
classifies  the  Q, ( i- 1 ,2 ,  •  •  •  ,n)  ,  and  computes  X3.  Estimators  of  the 
two  "true"  egression  lines  are  computed,  and  the  two  minimum  and  maximum 
sample  values  are  determined.  In  addition  to  the  printout  of  the  above 
computed  statistics,  optional  features  provide  for  printout  of  the 
categorization,  a  plot  of  the  data,  and  a  plot  of  the  data  with  contour 
ellipses . 
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II.  COMPUTATIONAL  PROCEDURE 


i he  estimators  of  the  five  pop  lation  parameters  are 
*  *v  =  ?v  =  sv  =  I(xvl  •^v)3/(n-l)  ,  (9) 


j  “Xj  )  (Xj  ,  -x^ ) 

VZ(XM-^>a  ZXXs,-^)2 


where  —1,2  and  all  summation  signs  run  from  1  to  n.  The  minimum 
and  maximum  sample  values  of  each  random  variable  are  the  respective 
order  statistics  Xj(r)  and  x^j),  Xj^j.  Estimates  of  the 

two  "true"  regression  lines  are 


Xj  =  X.,  +  r  (Xj-Xj)  =  a,  +  bjXj, 


x,  =  xt  +  r  J.  (xc-x-,)  =  a2  +  b^ 

s„ 


(10) 


The  five  statistics  of  equation  (9)  are  substituted  for  the 
five  parameters  of  equation  (2),  giving 


1 

*  2r  fxu-xi|  I 

|  *21  _X2  \  j.  1  *2  t  _X2^ 

2 

( l-r;  ) 

l  ! 

2r  !  1  ! 

i  S2  1  1  S2  1 

After  evaluating  §, (i=l ,2 ,  •  • • ,n) ,  the  classification  of  the  n 
observations  is  performed  by  comparing  each  Qt  with  (k-1)  specified 

c. -values  (see  card  type  5)  which  define  k  mutually  exclusive  and 

_  *  ^ 

exhaustive  areas  between  coaxial  contour  ellipses.  If  cz  s;  Q,  <  c‘ 

PJ-1  p3 

( j -l  ,2 ,  •  ■  •  ,k)  ,  the  i**1  observation  falls  within  the 
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t  h 

j  contour  ellipse  and  outside  the  (j-1)  contour  ellipse.  The 
number  of  (xt , ,Xg ,) -observations ,  0^ ( j=l,2, • *  * ,k) ,  falling  within 
each  of  the  k  areas  is  then  determined,  and  X2  as  given  by  equation 
(6)  is  computed. 

If  the  quantity  (k-6)  >  0,  the  statement  "CHI-SQUARE  = 

WITH  (k-6)  DEGREES  OF  FREEDOM"  is  printed  (see  page  19) ;  if  the 
quantity  (k-6)  s  0,  the  statement  "CHI-SQUARE  COULD  NOT  BE  CONFUTED 
BECAUSE  OF  INSUFFICIENT  NUMBER  OF  INTERVALS"  is  printed. 


III.  INPUT  PREPARATION 
A.  Deck  Setup 

The  input  deck  is  listed  below  by  card  type.  Multiple  jobs 
may  be  processed  by  stacking  sets  of  card  types  1  through  10.  Up  to 
twenty  runs  may  be  performed  for  each  set  of  sample  data  by  stacking 
sets  of  card  types  4  through  9  behind  card  type  10. 


CARD  TYPE  1  -  JOB  IDENTIFICATION  CARD 

CARD  TYPE  2  -  VARIABLE  FORMAT  CARD 

CARD  TYPE  3  -  NUMBER  OF  RUNS  CARD 

CARD  TYPE  4  -  MAIN  CONTROL  CARD 

CARD  TYPE  5  -  INTERVAL  IDENTIFICATION  CARD 

CARD  TYPE  6  -  CONTOUR  IDENTIFICATION  CARD 

CARD  TYPE  7  -  PLOT  HEADING  CARD  (optional) 

CARD  TYPE  8  -  TRANSFORMATION  CONSTANT  CARD 

CARD  TYPE  9  -  x2  TRANSFORMATION  CONSTANT  CARD 

CARD  TYPE  10-  SAMPLE  DATA  CARD 
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B. 

Input  Deck  Description 

CARD  TYPE  1 

-  JOB  IDENTIFICATION  CARD 

Column 

Format 

Program  Variable 

Explanation 

1-80 

10A8 

J0B(1)-J0B(1O) 

For  job  identification. 

CARD  TYPE  2 

-  VARIABLE  FORMAT  CARD 

Column 

Format 

Program  Variable 

Explanation 

1  80 

10A8 

fWT(l) -FMT(IO) 

The  format  by  which  the  (x^Xg)- 
values  (card  type  10  or  input 
record  if  data  is  on  tape)  are 
to  be  read.  The  format  specifi¬ 
cations  must  be  enclosed  by 
parentheses . 

CARD  TYPE  3 

-  NUMBER  OF  RUNS  CARD 

Column 

Format 

Program  Variable 

Explanation 

4-5 

12 

NRUN 

The  number  of  runs  to  be 
performed  per  set  of  sample 
data.  1  s  NRUN  s  20. 

CARD  TYPE  4 

-  MAIN  CONTROL  CARD 

Column 

Format 

Program  Variable 

Explanation 

5 

11 

NPRINT 

0  or  blank  -  Categorization 
is  not  to  be  printed. 

1  -  Categorization  ijs  to  be 
printed. 

10 

11 

JPL0T 

0  or  blank  -  No  plotting  is  to 
be  performed. 

1  -  Both  data  and  contours  are 
to  be  plotted. 

2  -  Data  only  is  to  be  plotted. 

14-15 

12 

N0TIC 

The  number  of  desired  tick  marks 

on  Che  axis  corresponding  to  the 
larger  range  of  xx  and  X£ .  If 
N0TIC  =  0  or  blank,  N0TIC  is  set 
equal  to  15. 
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Column 

Format 

Program  Variable 

Explanation 

20 

11 

INPUT 

0  or  blank  -  Sample  data  Is  on 
punched  cards . 

1  -  Sample  data  is  on  tape. 

25 

11 

NTYPES 

0  or  blank  -  Plot  headings 
(card  type  7)  are  not  to  be 
input . 

1  -  Plot  headings  are  to  be 
input . 

29-30 

12 

ITRAN1 

Number  corresponding  to  the 
transformation  to  be  performed 
on  the  xx -values.  1  £  ITRAN1  £  13. 

34-35 

12 

ITRAN2 

Number  corresponding  to  the 
transformation  to  be  performed 
on  the  Kg -values .  1  £  ITKAN2  £  13 . 

39-40 

12 

ISA 

The  number  of  decimal  digits 
to  be  used  in  labeling  abscissa 
tick  marks.  0  £  ISA  £  14. 

44-45 

12 

IR0 

The  number  of  decimal  digits 
to  be  used  in  labeling  ordinate 
tick  marks.  0  £  IR0  £  14. 

49-50 

12 

NPAIR 

The  number  of  pairs  of  (XjjXg)- 
values  per  sample  data  card  or 
sample  data  input  record  if 
data  is  on  tape. 

55 

11 

NN 

The  number  of  sets  of  probability 
values  to  be  used  to  define 
interval  bounds  for  the  Chi-Square 
Test.  1  £  NN  £  5. 

60 

11 

NL 

The  number  of  sets  of  probability 
values  to  be  used  to  define 
contours  to  be  plotted. 

1  £  NL  £  5. 
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TRANSFORMATIONS 


The  following  numbers  are  used  to  identify  the  thirteen  available 
transformations  (ITRAN1  and  ITRAN2) . 

Transformation  Number  Transformation 


1 

X  «■  X 

2 

x  4-  In  x 

3 

x  <-  ln(ln  x) 

4 

x  4-  ln(A+x) 

5 

x  4-  ln(B+ln(C+x)) 

6 

x  «-  /x 

7 

X  4-  1/x 

8 

x  4-  1/ (D+x) 

9 

x  4-  sin_1x 

10 

x  ♦-  2  sin  _1/x 

11 

x  «-  x/E 

12 

x  «-  sin  x 

13 

X  4-  COS  X 

Constants  A,B,C,D,  and  E  are  input  on  card  type  &  for  x1  and  card 
type  9  for  Xg . 
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CARD  TYPE  5 

-  INTERVAL  IDENTIFICATION  CARD 

Column 

Format 

Proarao  Variable 

Explanation 

1-5 

F5.3 

PS(1) 

The  smallest  probability  value 
for  the  1st  set  of  probability 
values  defining  interval  bounds 
for  the  Chi-Square  Test. 

6-10 

F5.3 

PDELT(l) 

The  increment  value  for  the 

1st  set  of  probability  values 
defining  Interval  bounds  for 
the  Chi-Square  Test. 

11-15 

F5.3 

PE(1) 

The  largest  probability  value  for 
the  1st  set  of  probability  values 
defining  interval  bounds  for  the 
Chi-Square  Test. 

16-20 

F5.3 

PS  (2) 

II 

21-25 

F5.3 

PDELT(2) 

II 

26-30 

F5.3 

PE  (2) 

II 

31-35 

F5.3 

PS  (3) 

It 

36-40 

F5.3 

PDELT(3) 

It 

41-45 

F5.3 

PE<3) 

II 

46-50 

F5.3 

PS  (4) 

II 

51-55 

F5.3 

PDELT(4) 

II 

56-60 

F5.3 

PE  (4) 

II 

61-65 

F5.3 

PS  (5) 

The  smallest,  increment,  and 

66-70 

F5.3 

PDELT(5) 

largest  probability  value  for  the 

71-75 

F5.3 

PE  (5) 

5th  set  of  probability  values 
defining  interval  bounds  for  the 
Chi-Square  Test. 

CARD  TYPE 

6  -  CONTOUR  IDENTIFICATION  CARD 

This 

card  type  is 

omitted  if  JPL0T  = 

0  or  JPL0T  *  2  on  card  type  4. 

Column 

Format 

Proaram  Variable 

Explanation 

1-5 

F5.3 

PPS(l) 

The  smallest  probability  value 

for  Che  1st  set  of  probability 
values  defining  contours  to  be 
plotted . 
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Column 

Format 

Program  Variable 

Explanation 

6-10 

F5.3 

PPDELT(l) 

The  increment  value  for  the 

1st  set  of  probability  values 
defining  contours  to  be  plotted. 

11-15 

F5.3 

PPE(l) 

The  largest  probability  value 
for  the  1st  set  of  probability 
values  defining  contours  to  be 
plotted. 

16-20 

F5.3 

PPS(2) 

ft 

21-25 

F5.3 

PPDELT(2) 

II 

26-30 

F5.3 

PPE(2) 

II 

31-35 

F5.3 

PPS(3) 

It 

36-40 

F5.3 

PPDELT(3) 

II 

41-45 

F5.3 

PPE(3) 

II 

46-50 

F5.3 

PPS(4) 

II 

51-55 

F5.3 

PPDELT(4) 

II 

56-60 

F5.3 

PPE(4) 

fl 

61-65 

F5.3 

PPS(5) 

The  smallest,  increment,  and 

66-70 

F5.3 

PPDELT(5) 

largest  probability  value  for 

71-75 

F5.3 

PPE(5) 

the  5th  set  of  probability 

values  defining  contours  to  be 
plotted. 


CARD  TYPE  7  -  PLOT  MAD  IMG  CARD 

This  card  type  is  omitted  if  NTYFES  ■  0  on  card  type  4.  Two  card 
types  7  are  used  for  identification. 


Column 

Format 

Program  Variable 

Explanation 

1-80 

1QA8 

XID1(1)  -  XID1(10) 

The  first  line  of 
identification. 

the  plot 

1-80 

1QA8 

XID2(  1)  -  X ID 2 (10) 

The  second  line  of 
Identification. 

the  plot 

This  card  type  is  omitted  if  ITRAN1  -  l  on  card  type  4.  If 
2  ••  ITRAN1  v  13,  but  the  transformation  corresponding  to  the  transformation 
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number  does  not  contain  one  of  the  constants  A ,8,0,0,  or  E,  a  blank  card 
must  replace  card  type  8. 


Column 

Format 

Program  Variable 

Explanation 

1-14 

E14.8 

ATRAN1 

The  constant 
number  4  for 

A  in  transformation 
transforming  xx -values 

15-28 

E14.8 

BTRAN1 

The  constant 
number  5  for 

B  in  transformation 
transforming  xx -values 

29-42 

E14.8 

CTRAN1 

The  constant 
number  5  for 

C  in  transformation 
transforming  Xj -values 

43-56 

E14.8 

DTRAN1 

The  constant 
number  8  for 

D  in  transformation 
transforming  Xj -values 

57-70 

E14.8 

ETRAN1 

The  constant 

E  in  transformation 

number  11  for  transforming  -values . 

CARD  TYPE  9  -  *a  TRANSFORMATION  CONSTANT  CARD 
This  card  type  is  omitted  if  ITRAN2  =  1  on  card  type  4.  If  2  s  ITRAN2  s  13, 
but  the  transformation  corresponding  to  the  transformation  number  does  not 
contain  one  of  the  constants  A,B,C,D,  or  E,  a  blank  card  must  replace  card  type  9. 


Column 

Format 

ProRram  Variable 

Explanation 

1-14 

E14.8 

ATRAN2 

The  constant 
number  4  for 

A  in  transformation 
transforming  -values 

15-28 

E14.8 

BTRAN2 

The  constant 
number  5  for 

B  in  transformation 
transforming  x^ -values 

29-42 

E14.8 

C TRAN 2 

The  constant 
number  5  for 

C  in  transformation 
transforming  Xg -values 

43-56 

E14.8 

DTRAN2 

The  constant 
number  8  for 

D  in  transformation 
transforming  x^ -values 

57-70 

E  14.8 

ETRAN2 

The  constant  E  in  transformation 
number  11  for  transforming  x,- 

values . 
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CARD  TYPE  10  -  SAMPLE  DATA  CARD 


Column 

Format 

Program  Variable 

Explanation 

Variable 

Specified 
on  card 
type  2 

IEND 

Indicates  termination  of  the 
sample  data  cards.* 

Variable 

Specified 
on  card 
type  2 

X1(I)  &  X2(I) 

The  ith  pair  of  (x^.Xg,)- 
values;  1  s  4,000. 

*If  IEND  * 

0  or  blank, 

NPAIR  (xx .Xg) -values 

are  on  the  card  and  at  least 

one  more  sample  data  card  follows. 

If  IEND  >  0,  IEND  (xx .Xg) -values  are  on  the  card  and  this  card  is  the 
last  sample  data  card. 


C .  Request  Sheets 


If  the  sample  data  (card  type  10)  is  input  on  punched  cards , 
the  job  request  sheet  is  prepared  as  shown  below. 


SffURITY  CLASS  1  T  l(  ATI  ON 

□  ts  r~is-  rir.  riiu. 

N  AMf 

J.  BROWN 

FBI 

llffC.  ~ 

1200 

BSEjam 

CBS 

BiH 

[  1  rnwp | L  t  (X  1  GO 

9 

1 _ IQB  CARD _ 1 

mma 

XI 

(MARC 

it  CODE  i 

snrnm 

1  BI  -  CHI 

1  |  COM* 11  GO 

u 

5 

8 

B 

c 

A  3 

“1 

E9HH 

TAPES  CALLED 

FOR  BY  PROBLEM  PROGRAM 

TAPI  Nl'MPf  P 

as 

h 

- , 

FI  LC  p  *<'  Tt  r  T  *iN 

RROCiftAMMl  N  Nl'WBI  R 

j 

SP  (  C 1 *1  MANPL  1 NG 

*tt*ch*d  gnat.) 

j 

0»fRAT0R*S  rcuMfNTS 


[J...OJ  □hold  Qi 


□  IT  CHlCKT*.  ill  RC  Vf  RSI 
Slot  Tor  ftoninr-NAi 


S*  f  C  I  ft  L  INS  TRLir  T  IONS  (  con  T  i  nurd  on  rrvrrtr) 


7010  JOB  REQUIEST  NDW.NWL  -1210  i*  IRCV.  tl.AAl 


If  the  sample  data  Is  inpuc  on  tape,  the  tape  number  is  written  on 
the  Job  request  sheet  instead  of  "Scratch"  and  the  tape  number  is 
punched  in  place  of  XXXX  on  the  "UEL,  PUUOOOt"  card  (second  card) 
of  the  IOD  deck. 
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OATE  AND  TIME  PROBLEM  RAN 


16 


TOO 

in 


II 

cvl 

X 


Su  S  aS  m  X  X 


17 


CHI-SQUARE  TEST 


M 


ot 

G 

W-t  CV 

<Cd 

C9 

41 

o 

<  W  <Bd 

<w 

1*4 

<w 

o 

\ 

| 

Cl  0) 

Cl 

M 

Cl 

S3 

<¥ 

#•  ••  rH 

| 

% 

1 

CV 

<w  <c y 

« •  ••  cgp 

M 

<f9 

r-i  CV 

X 

M 

53 

o  o 

o 

O 

o 

O 

w  w 

w 
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B.  CRT  Output 


The  plots  are  performed  by  the  system  "GRP  Plot  Subroutine"  and 
are  plotted  by  the  CRT  printer  units.  The  plots  are  in  an  ordinary 
Cartesian  coordinate  system,  and  the  origin  is  at  (x^j  • 

The  abscissa  of  the  system  is  the  x^axis;  the  ordinate  of  the  system 
is  the  Xg-axis.  An  asterisk  (*)  is  used  to  denote  an  (x1( ,*3,)- 
observation;  a  zero  (0)  is  used  to  denote  the  center  (xj,Xg)  of  the 
distribution  of  the  observations.  In  addition,  a  zero  (0)  is  used 
to  denote  points  on  each  of  the  estimated  regression  lines, 
x^  =  f(xx)  and  xx  =  f(xa) .  Up  to  four  points  (depending  upon  the  range 
of  the  particular  plot)  are  plotted  for  each  regression  line.  The 
plotted  points  on  the  regression  lines,  however,  are  not  connected. 
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C .  Optional  Output 


The  output  illustrated  on  pages  IS,  16,  17,  and  19  of  Section  IV .A 
is  automatically  printed  for  each  run  processed.  The  categorization 
(page  18)  and  the  CRT  output  is  controlled  by  the  analyst.  This 
flexibility  is  provided  by  the  program  variables,  NPRINT  and  JPL0T, 
on  card  type  4.  The  allowable  input  values  (two  for  NPRINT  and  three 
for  JPL0T)  for  these  program  variables  permit  six  different  output 
options.  A  quick  reference  of  the  available  output  options  is  given 
in  the  following  table. 


NPRINT 

JPL0T 

OUTPUT 

0 

0 

1. 

Chi-Square  Test 

0 

1 

1. 

Chi-Square  Test 

2. 

Plot  of  data  with  contours 

0 

2 

1. 

Chi-Square  Test 

2. 

Plot  of  data 

1 

0 

1. 

Chi-Square  Test 

2. 

Categorization 

1 

1 

1. 

Chi-Square  Test 

2. 

Categorization 

3. 

Plot  of  data  with  contours 

1 

2 

1. 

Chi-Square  Test 

2. 

Categorization 

3. 

Plot  of  data 
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.  Program  Running  Time 


I 


The  running  time  is  printed  at  the  end  of  each  run  processed. 
The  following  table  of  observed  times  may  be  used  as  a  guide  for 
estimating  program  running  times. 


Sample  Size 

No.  of  Intervals 

Categorization 

Plots  Time 

in  Seconds 

50 

10 

Yes 

Data  & 
Contours 

8.2 

100 

14 

No 

Data  & 
Contours 

9.8 

150 

12 

Yes 

Data  & 
Contours 

14.3 

500 

100 

No 

No 

16.2 

1000 

150 

Yes 

Data  & 
Contours 

56.5 

2500 

200 

No 

Data  & 
Contours 

81.2 
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V.  NUMERICAL  EXAMPLE 


A.  Problem  Description 

The  example  problem  consists  of  100  observations  (pairs  of  (Xj.Xg)- 
values)  randomly  sampled  from  a  continuous  two-dimensional  distribution. 

Two  runs  were  performed,  one  on  the  original  observations  and  one  on  the 
logarithmically  transformed  values.  Because  some  of  the  original 
observations  (both  x^ -values  and  x^ -values)  were  negative,  transformation 
number  4  was  used.  For  transforming  x1 ,  the  constant  A  was  set  equal  to 
25;  for  transforming  Xg ,  A  was  set  equal  to  50.  In  both  runs,  14  intervals 
were  specified  for  the  Chi-Square  Test.  Categorisation  was  requested  for 
the  first  run,  but  it  was  not  requested  for  the  second  tun.  The  heading 
of  the  plot  (card  type  7)  was  input  for  the  first  run  only. 

The  input  deck  setup  is  illustrated  on  the  following  Data  Card 
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CHI—  SQUARE  TEST 

*U*  NUMBER  1  WITH  100  OBSERVATIONS  AND  14  INTERVALS 
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D.  Discussion  of  Test  Results 


Suppose  the  0.05-level  of  significance  is  the  preselected  level  for 
testing  the  null  hypothesis  of  normality.  To  test  the  null  hypothesis 
"The  100  randomly  sampled  observations  are  from  a  parent  bivariate 
normal  population,"  Xs  =  21.5  (from  the  first  run)  is  compared  with 
y2(0.95,8)  =  15.5.  The  null  hypothesis  is,  therefore,  rejected  at  the 
0.05-level  of  significance.  To  test  the  null  hypothesis  of  normality 
of  the  transformed  random  variables,  Xs  =  13.5  (from  the  second  run)  is 
compared  with  ^(O^.S)  =  15.5.  Because  the  latter  null  hypothesis 
cannot  be  rejected  at  the  0.05-level  of  significance,  we  conclude  that  the 
two  original  random  variables  are  approximately  jointly  log  normally 


distributed. 


VI.  FLOW  CHART 
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VII.  PROGRAM  LISTING 


DIMENSION  TABL.EC  ( 999  )  «  JOB  (  1 0  )  <  0BS1999) ,XlORI6(4000 ) ♦ 

|  X20R 16(4000 ) »PVAL (10) .PPVAL ( 1 0 1 »TABP ( 10 ) *TABN( 10) 

COMMON  P  S  (1 0  )  , POELT (10).PE(10).PPS(10). PPOEL  T ( 1 0  > « PPE ( 1 0  >  » 

1  XI (4000 ) «X2(4000 ) «P (999 ) »C2 (999 ) *0(4000 ) ♦ HANDY (  18 )  * 

2  FMT ( 1 0 ) «X 101 (10) .  X I D2  (  1  0  )  «  1 OBS (999 ) «PP (999 ) *E (*000 ) .C2P (999 ) 
DATA  (TABP( I),I"1.10)(1..10.«100. .1000. . 10000* « 1 00000. . 1 000000 •» 


1  10000666. 
DATA (TABN( I ) . I ■ I , 10) (.606006001 

1  .00001 . .0001 ..001 . .01 

DATA  (TABLECd  ).  1-1.36) 

1 (.20010007E-02.  .40040053E-02, 

2  . 1002S0B4E-01 .  . 12036 145E-01 . 

3  • 18081 489E-01 .  .20100672E-01 

4  .26 1 70479E-0 1 .  .281 97849E-0 1 

5  .3429231  BE -01.  .36327941E-01 

6  .42447273E-01 .  .44491218E-01 

7  .5063561 6E-0 1 .  .52687951 E-0 1 

8  • 58057621 E-0 1.  .6091 84 1 5E-0 1 

9  .671 1 3567E-01 .  . 691 82890E-0 1 
DATA  (TABLECd  ).l  >37 .72) 

1 (.75403734E-01 .  .77481657E-01 

2  .83728408E-01 .  .8581 5002E-0 1 

3  .92087877E— 01 .  • 94 1 832 1 SE-0 1 

4  • 10048243E— 00.  • 102586S9E-60 

5  . 1 089 1 237E-00 .  .11102542E-00 

6  .1 1737799E-00.  • 1 195000 1 E-00 

7  • 1 2S87960E— 00 .  . 12801 066E-00 

8  . 13441 750E-00,  • 1 3655768E-00 

9  • 1 4299200E— 00 .  . 14514139E-00, 
DATA  (TABLECd  ).  I  -73. 108) 

1 (. 15160343E-00.  • 1 537S209E-00 

2  .160252096-00.  . 1 62420 1 1 E-00 

3  .16893831 E-00.  , 1 71 1 1 578E-00 

4  . 1 7766243E-00 .  . 1 7984942E-00 

5  • 1 B642476E— 00 .  . 1 8862 1 36E-00 

6  • 1 9522566E -00 .  . 1 9743 1 95E-00 

7  . 20406545E-00,  .206281 52E-00 

8  .21 294449E— 00 .  .2 1 51 7042E-00 

9  .2218631 2E-00.  .22*0990 1 E-00 
DATA  (TABLECd  ).  I  "109.  144) 

1  (.23082170E-00.  . 23306763E-00 

2  .23982059E-00.  .24207666E-00 

3  •  248860 16E-00.  .25U2645E-00 

4  • 257940 76E— 00 .  .26021 737E-00 

5  .26706279E-00.  . 2693498 1 E-00 

6  .27622660E-00.  . 278524 l 3E-00 

7  .28543260E-00.  .28774074E-00 

8  .294681 18E-00.  . 29700002E-00 

9  .3039727 IE-00.  . 30630236E-00 
DATA  (TABLECd  ).  1-145,  180) 

1 t.3l330762E-00,  . 3 » 5648 1 7E-00 

2  . 322686 3 OE -00 «  . 32503786E-00 

3  .3321 091 7E-00.  .33447164E-00 
A  . 34  I S766AE -00  »  .34 395053E-00 
:  .351 0891 5E-00,  . 3534 7436E-C 0 

6  .3606471  IE-00,  . 3630«375E-00 

7  . 37025097E -00 .  . 372659 1 6E-00 

8  .379901 1 7E-00,  . 38232 1 0 1 E-00 

9  . 389598 16E -00.  . 39202977E-00 


1 06000006. . 1 600066000.  ) 

•  00000001 1 .0000061  * .600001 ♦ 
•1.1.) 

60090 180E-02.  .80160428fc-02« 
14049230E-01 *  .160643436-01, 
22121895E-01 .  .24145162E-0 1 . 
30227276E-01 .  .32258764E-0 1 . 
38365639E-01 .  .40*0541 5E-0 1 « 
46537254E-0I «  .48585385E-0 1 • 
54742394E— 01 .  .56798949E-0 1 . 
6298 1 334E— 0 1 .  .65046383E-0 1 . 
71254355E— 01  *  .73327969E-0 1  ) 

79561740E— 01 .  .81643989E-0 1 . 
87903775E— 01 .  . 89994 732E-0 1 . 
96280751E— 01 .  .98380488E-0 1 . 
1 0469296E— 00.  . 10680 155E-00, 
11314070E-00.  . 1 1525823E-00. 
121 6242BE— 00  «  . 1 237508 1 E-00 . 
1301 4399E— 00  »  . 1 322796 1 E-00 . 
13870016E— 00.  . 14084493E-00. 
1 4729308E— 00 .  . 14944709E-00 ) 

1 5592308E— 00  «  . 1580864 1 E-00. 
16459049E— 00.  . 1 6676322E-00 , 
17329561E-00.  . 17547783E-00. 
1 8203880E— 00 •  . 18423058E-00. 
19062037E-00.  . 19302 1 80E-00. 
19964067E— 00*  .20185184E-00. 
20850004E— 00,  .21072103E-00. 
217398B3E-00.  .21962973E-00. 
22633740E-00.  .22B57829E-00 ) 

1 2353 1 609E— 00 .  .23756707E-00, 
24433527E— 00.  , 2465964 3E-00 . 

1 2533953 1 E-00 .  .255666746-00. 
1 26249657E-00.  .26477838E-00. 

1 27  I 63945E-00.  .273931 7 1 E-00 . 
1 28082431 E-00.  .2831271 3E-00, 
,290051546-00.  . 29236502E-00 « 
1 29932 1 55E-00 .  .301645786-00. 
1 30863472E— 00.  . 3 1 09698 1 E-00 ) 

1 3 1 799 1 46E— 00.  .32033750E-00. 
i  3273921 9E-00.  . 32974929E-00, 

>  3 368 3 73 OE- 00 .  . 33920557E-00 , 

>  34632724E -00 •  .34870677E-00. 
>355862*2E-00«  . 3S825333E-00 . 
>3654432 7E -00.  . 36784568E-00 . 
,375070256-00.  .377484256-00. 
,  38474379E— 00 .  . 387 1 6950E-00 • 

>  394 4 64 34 E-00 .  . 39690 1 88E-00 ) 


8C000010 
BC000020 
BC000030 
BC000040 
BC 000050 
BC000060 
BC000070 
BC000080 
BC000090 
BC000100 
BCOOOl 10 
BC000120 
BCOOOl 30 
BC000140 
BCOOOl 50 
BC000160 
BCOOOl 70 
BCOOOISO 
BCOOOl 90 
BC000200 
BC000210 
BC000220 
BC000230 
BC000240 
BC000250 
BC000260 
BC000270 
BC000280 
BC000290 
BC000300 

ecooosio 
BC000320 
BC000330 
BC000340 
BC000350 
BC000360 
BC000370 
BC000380 
BC000390 
BC000400 
BC000*»10 
BC000420 
BC000430 
BC000440 
BC000450 
8C000460 
BC000470 
BC000*80 
BC 000*90 
BC000500 
BC000510 
BCC00520 
BC000530 
BC000540 
BC000550 
BC000560 
HC000570 
BC0005H0 
HC000590 
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DATA  (TABLEClI) 

1 < . 399342398-00  « 

2  *4091 3433E-00, 

3  «4 1 897445E -00 • 

4  .42886322E-00, 

5  .438801 13E-00. 

6  .448788678-00. 

7  .45882633E-'  0 . 

8  .4689 1 462E -00 . 

9  •  47905406E -00 . 
DATA  (TABLEC(l) 

1  (  .489245 1 7E-00  . 

2  . 49948847E -00 . 

3  . 50978450E+00 . 

4  .520133818+00, 

5  *  53053696E+00 , 

6  54099450E+00. 

7  .551507008+00. 

8  . 56207506E+00  , 

9  .572699258+00, 
DATA  (T1Bl.EC  11) 

I  . 583380 19E+00. 
_  .5941 1847E+00. 

3  .6049 1 472E+00 . 

4  .61 576956E+00. 

5  .626683648+00. 

6  .63  765760E+00 « 

7  .6486921 1E+00. 

8  .659787848+00, 

9  .670945478+00, 
DATA  (TABLECU) 

1  (  .682 16570E  +  0G, 

2  .693449238+00. 

3  .704796778+00. 

4  .716209078+00, 

5  . 72 76868 7E  + 00  » 

6  . 7392309 1E+00, 

7  ,750841978+00. 

8  . 76252084E  +  00 « 

9  .774268308+00, 
DATA  (TABtEC(I) 

1  (.786085188  +  00, 

2  .79797228E+00, 

3  • 80993G478+00 , 

4  .821960588+00, 

5  .83406349E+00. 

6  • 8462« 009E +00 , 

7  .858491278+00, 

8  .870817978+00. 

9  .883221118+00, 
DATA  (TABlEC(I) 

1 < .895701658+00, 

2  • 90826056E+C0 , 

3  .720898838+00, 

4  ,9336 174QE+00, 

5  .946417528+00, 

6  .959300018+00. 

7  ,  97226602E  +  00 , 
.8  .985316648+00. 


1-181 ,216) 

.401 78588E-00 

.411589838-00 

.421442068-00 

.431 34307E-00 

.441293348-00 

•45129336E-00 

.461343648-00 

.471444678-00 

.481596978-00 

1-217. 252) 

•491801 08E-00 

.502057518+00 

.512366818+00 

.522729538+00 

•  533 1 46228+00 

.543617458+00 

. 554 1 4379E+00 

. 56472582 E+ 00 

.575364148+00 

1-253,288) 

.586059368+00 

.596812078+00 

.607622918+00 

.618492508+00 

. 62942 14 9E+ 00 

.640410538+00 

.651460288+00 

. 66257 142E+00 

. 67374463E+00 

I -289,324  ) 

.68498062E+00 

. 69628008E+00 

.707643758+00 

. 71907235E+00 

. 73056664E+00 

.742127368+00 

,753755308+00 

.765451248+00 

.777215988+00 

1 -325,360) 

.789050348+00 

. 800955 1 3E+00 

.812931228+00 

.824979458+00 

. 8371007 0E +00 

. 84929586E+00 

.861565838+00 

.87391 155E+00 

.886333958+00 

1 >361 ,396) 

.898833998+00 

.91 1412658+00 

.924070928+00 

.936809828+00 

.949630378+00 

.962533648+00 

.975520708+00 

.988592648+00 


40423237E-00 
414048348-00 
42391 272E-00 
43382600E-00 
44378866E-00 
453801 20E-00 
463864118-00 
47397792E-00 
4841 431 2E-00 

49436026E-00 
50462986E+00 
51495246E+00 
52532862E  +00 
53575889E+00 
546243B4E+00 
55678405E+00 
5673801 OE+OO 
57803259E+00 

5887421 2E+00 
59950931 E+00 
61033477E+00 
621 2 1 915E+00 
632 1 6309E+00 
6431 6725E+00 
65423228E+00 
6653588BE+00 
67654772E  +00 

68779950E+00 
6991 1495E+00 
71049478E+00 
72193974E+00 
733450S6E+00 
74502802E+00 
756672B6E+00 
76838595E+00 
780168018+00 

79201990E+00 
803942448+00 
8159364BE+00 
82800288E+00 
840142528+00 
8 52 3563 OE+OO 
064045 1 2E+00 
87700992E+00 
889451 64E+00 

90197125E+00 
914569718+00 
92 72 4 804 E+00 
94000726E  +00 
952848398+00 
96577251E+00 
97878069E+00 
991 87402E  +00 


406681858-00. 
4 1 6509888-00 « 
42638644E-00 , 
436312028-00. 
446287108-00. 
45631219E-00, 
466387778-00. 
476514388-00. 
486692528-00 ) 

49692272E-00. 
507205528+00. 
5 1 754 1 46E+00 , 
527931098+00. 
538374988+00. 
548873698+00 . 
55942781 E+00. 
57003791 E+00. 
58070460E+00 ) 

591428498+00. 
<60221 019E+00. 
61 305032E+00, 
6239495 3E+00. 
6349084 6E+00, 
64592777E+00, 
657008 13E+00. 
668150228+00, 
67935474E+00 > 

690622378+00. 
701953858+00, 
713349898+00. 
72481 1248+00. 
736338658+00. 
747932888+00. 
759594728+00. 
771324968+00. 
7831244  18+00  ) 

794993888+00, 
80693421 E+00, 
018946268+00. 
831030898+00, 
H4318t'98e+00. 
055421 431+00. 
8677291 78+00. 
0801131 18+00. 
892574218+001 

9051 13438+00, 
917731778+00. 
930430238+00. 
943209S2E+00 . 
956071608+00. 
969016638+00. 
982045998+00. 
995160798+00. 


BC000600 

BC00061 0 
BC000620 
BC000630 
BC000640 
BC000650 
BC000660 
BC000670 
BC000680 
BC000690 
BC000700 
BC000710 
BC000720 
BC000730 
BC000740 
BC000750 
BC000760 
BC000770 
BC000780 
BC 000790 
BC000800 
BC000810 
BC000820 
BC000830 
BC000840 
BC000850 
BC000860 
BC000870 
BC000880 
BC000890 
BC000900 
BC000910 
BC000920 
BC000930 
BC000940 
BC000950 
BC000960 
BC000970 
BC000980 
BC000990 
BC001000 
BC001010 
«C0010?0 
BCOOI 030 
BC001040 
UC001050 
BC001060 
BC001070 

bcooiobo 

8C001090 
8C00U00 
8C001 I  10 
BC001 120 
BCO01 130 
BC001 140 
BC001 ISO 
BCOOI 160 
BCOOl 170 
BCOOI 180 
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9  .99845298E+00, 
DATA  ( TABlEC  <  I  ) 
I (.101 1 6762E+0 1 . 

2  . 10249874E  +  01  , 

3  *  1 0383877E+0 1  . 

4  . 1051878SE+01 . 
b  . 10654609E+01  . 

6  . 10791 362E+01  . 

7  . I0929056E  +  01  . 
6  « 1 1067705E+0I . 
9  • 1  1207321E+01  . 

DATA  (TABLEC(I) 

1  (.1 1347920E  +  01 « 

2  • 1 148951 3E  +  01  • 

3  . 1 16321 1 6E  +  0 I « 

4  , 1 1 775743E+01 . 

5  .1 1 920409E+0 I , 

6  •  12066130E+01 • 

7  .  122129I9E+01 • 

8  • 12360794E+01  . 

9  . 12509771E+01 . 
DATA  (TABlEC(I) 

1 (. 12659865E+01 • 

2  . 1281 I095E+01 . 

3  •  12963476E+01 . 

4  . 131 17028E+01 . 

5  . 13271 768E+01  . 

6  . 13427714E+01  . 

7  .  I3584886E  +  01  . 

8  . 13743302E  +  01  . 

9  . 1 3902984E+0 1 « 
DATA  <TABcEC(ll 

1  (.14063950E+01  « 

2  .  M226223E  +  01 . 

3  .  1  4389823E  +  0  1  » 

4  .  14554773E  +  01  . 

5  .  14721094E.01  , 

6  .  14888809E+01  . 

7  . 15057944E+01 . 

8  . 15228520E+0I  . 

9  ,15400564001  , 
DATA  (TABLEC(l) 

1  (. 1 5574 l  0 l E  +  0 1 • 

2  . 15749 157E+01 . 

3  . 15925759E+01 . 

4  . 16103934E+01. 

5  • 162837 10E+0I . 

6  .164651 1 7E  +  01 , 

7  . 16648 185E *01 , 

8  . I6832944E+01 . 

9  . 17019425E+01 . 
DATA  ( T  ABLEC ( I  ) 

1  (. 17207662E+0I • 

2  .  I 7397607E  +  OI  , 

3  . 1 7589535E  +0 1 . 

4  •  1  778324 1 E  +  0 1 • 

5  .  17978842E+01  . 

6  •  181 76374E+01  , 

7  .  18375877001  , 


.  1 00 1 7506E+0 1 . 
. 1 *397.432  ) 

.  10149957E  +  OU 
. 1 028329 1 E  +  0 i . 

•  1 04 1 75 1 9E  +  0 1 . 
. I0552655E+01 . 
. 10688710E+01  . 
. 1 0825697E+0 1  . 
. 1 0963628E+0 1  . 
. 1 1 102518E+01 , 
. 1 1 242378E+0  1  * 

. I *433.468  > 

. 1 1 383224E+0 1 . 
. 1 1525C69E+01 . 
. 1 1 667926E+0 1 . 
, 1 181 1 8 1 2E  +  0 !  . 

•  1 1956740E+01 « 
.  12102726E  +  ''  1  , 
. 12249786E+  Jl , 
.  12397934E  01 . 
.125471 39E  01  » 

. I *469.504 
.  12697565c. +0 '  . 
. 1 2849081 E+0  . 
. 1 300 1 754E+0  * » 
. 1 3 1 5560 1 E+0 1  . 
.  1331064 OE+01  , 
.  1346689 1 E+0 1  . 
.  1  3624372E+0 1 « 
. 13783103E+01. 
. 1 39431 04E+0 1 < 
. 1*505.540) 

. 14104395E+0I  . 
. 14266998E+01. 
.  14430933E+01 . 
. 14596223E+01 . 
. 14762891 E+0 1 . 
. 14930959E+01 . 
. 15100452E+01 . 
,15271 393E+01 « 
. 1 5443808E+0 1 » 
•  I *541 ,576) 

, 1 561 7722E+0 1 • 
, 1 57931 62E  +  0 1  * 
. 15970154E+01 . 
. 16148727E+0! . 
. 16328908E+01 . 
. I  65 1 0727E  +  Q 1  , 
. 1 66942 1 5E  +  0  I  , 
. 1 6879401 E  +  0  1  « 
.  1  706631 9E  +  0 1 . 
. 1 *577,612  ) 

.  1 7254999E+0 1  * 

•  1 7445477E  +  0  1  , 
.  1 7637786E+0 1  . 
. 1 7831 962E+0 1 , 
. 1 8028042E  +  0  1  « 
.  1 8226064E  +  0  l  « 
. 1 8426065E  +  0  !  . 


.  1  0050536E  +  01  , 

.  10183207E+01 
. 1 031 6763E+01 

. 1 045 121 8E+0 1 
. 1 0586582E+0 1 
. 10722869E+01 
. 10860090E+01 

•  10998260E+01 
.  1 1 137391E+01 

•  1 1277497E  +  01 

.  1  1410591E+O1 
. 1 1560687E+01 
.1 1703801E+01 
. 1 1847946E+01 
.1 1993137E+01 
. 12139390E+01 
. 12286720E+01 
. 12435144E+0! 

. 12584677E+01 

. 12735337E+01 
. 12887140E+01 
. 13040105E+01 
• 13194248E+01 
. I3349589E+01 
. 13506145E+01 
. 13663937E+01 
. 13822984E+01 
. 13983305E+01 

• 14144922E+01 
■  1 4307856E+0 1 
. 14472128E+01 
.14637760E+01 
. 14804776E+01 
. 14973 198E +01 
. 15143050E+01 
, 1 531 4357C+0 1 
. 1 54871 45E+0 I 

• 15661438E+01 
. 15837263E+01 
. 1 601 4648E+01 
. 1 6 l 93620E+0 1 
. 1637420. E+01 
. 16556442E+01 
. 16740351E+01 
. 1692596 7E+01 
. 171 13322E+01 

. 17302449E+01 
. 17493381E+01 
. 1 7686154E+01 
. 17880802E+01 
. 18077364E+01 
. 18275877E+01 
•18476 380E +01 


. 1 00836225  +  0  1  ) 

.1021651 2E+0 1 « 
. 1 0350292E+0 1 » 

. 10484973E+0 ' • 

•  1 0620567E  +  0 1 , 
. 1 0757086E+0 1 . 
. 1 0894544E+  01. 
.11 032952E+0 I » 
.11! 72326E+0 1 , 
. 1 1312677E. 31 ) 

. 1 1454021E+01 . 
.1 1596370E+01 
.1173974 0E+01, 
. 1 1B84145E+01 . 
• 12029600E  01, 

•  1 21 761 2 1E+0 1 « 
. 1 232372 3E+0 1 . 
• 12472422E+0 1 . 
. 12622236E  +  01  i 

. 12773180E+01 
. 12925272E+0 1 « 
. I 3C78529E+0 1 . 

•  1 3232970E  +  0  1  « 
. 1 338861 3E +01. 
. 1 3545477E+0 1 » 
. 1 3703580E+  0 1 « 

•  1 3862944E+0 1 « 

. 1 4023587E+0 1  ) 

.141 8552  E+01, 
. 1 434B797E+0 1 • 
. 1 45 1 3407E+0  1  « 
, 1 4679384E+0 1 . 
. 14846748E+0  1  . 

•  1 50 1 5526E+0  1 . 
.151 85740E  +  0  1  . 
. 15357415E+01 . 
. 1 55305762  +  0  1  ) 

. 15705249E  +  01  . 
■  1 588 1 462E+0  1 , 
.1605924  1 E  +  01  , 

•  1 623861 4E  +  0  I  . 
.1641961 1 E  +  0 1  , 

•  1660226IE  +  0  1  « 

•  16786594E  +  01  . 
. 1 697264 2E  +  0  1  , 
. 17160436E+0  1  ) 

.1735001 1E  +  01  . 

•  1 754 1 400E  +  0 l . 

•  1 7734639E  +0  l  ♦ 
. 17929762E+01 . 
.181268086  +  01  . 
. 183258 15E+01. 
. 1652682 1E  +  0  I  . 


fctCOOl 190 
BC001200 
6C001210 
BC001220 

BC001230 
3C001240 
80001250 
BC001260 
BC001270 
BCO"'  1280 
dCO  4  290 
acoo i 3oo 
BC001310 
BC001320 
BC001330 
PC001340 
B  -.00 1350 
B  001350 
Bi_00  1  3-*0 
BC0013  0 
BCOO 1390 
B  001400 
B^OO 1410 
BC001420 
t  C00143C 
BC001440 
BCOOl  >50 
BC001460 
BCOOl 4 70 
BCOOl 480 
BC001490 
B '00 1 500 
BC001510 
BCOOl 520 
BCOOl 530 
BCOOl 540 
BCOOl 550 
BCOOl 560 
BCOOl 570 
BCOOl 580 
BCOOl 590 
BCOOl 600 
BC001610 
BCOOl 620 
BC00163C 
BCOO 1640 
BC001650 
BCOO 1 660 
8C00 I  670 
BCOOl 680 
BCOO 1 690 
BC001700 
BCOOl 710 
BCOO 1 720 
BCOOl 730 
BCOO! 740 
BCOOl 750 
BCOOl 760 
BCOOl 770 
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8  .18;  '7390E  +  0 1  ♦ 

9  . 18780954E+01 • 
DATA  < TABLEC ( I  ) 

1 (. 189866 12E+01 . 

2  .191 94406E+0 1 • 

3  • 19404381E+01 • 

4  . 1 96 1 6585E+0 1 « 

5  . 19831 064E+01  . 

6  .20047869E+01 . 

7  .20267049E+01 . 

8  .20488658E+01  . 

9  .2071 2750E+C 1  « 
DATA  < T  ABLEC  (  I  ) 

1 ( .2093938 1  E  +  0 1 . 

2  .21 16861 0E+01 • 

3  .21400497E+01 . 

4  .21 6351 03E+0 1 . 

5  .21B72495E+01 . 

6  .221 12738E+01 . 

7  .22355902E+01 « 

8  .22602059E+0 1 • 

9  .22851 284E+01 . 
DATA  (TABLEC(l) 

1 ( .231 03653E+0 1 . 

2  .23359247E+01 • 

3  .2361 8 1 S1E+01 « 

4  .23880449E+01  . 

5  . 24 1 46234E+0 1  * 

6  .244155986+01 . 

7  .24688640E+01 « 

8  .2496546 1 E+0 1 « 

9  .25246 1 68E  +  0 1  * 
DATA  ( TABLEC (I) 

1 ( . 25530870E+0 1 . 

2  .258 1 9684E+0 1 « 

3  .261 12729E+01 . 

4  .264 1 0 1 32E+0 1 . 

5  .2671 2025E+0 1 « 

6  .2701 8544E+0 1 • 

7  .27329825E+0 1 • 

8  • 27646047E+0 1 . 

9  .27967339E+01 , 
DATA  1 TABLEC (I) 


.  18628087E  +  01 
. 18832171E+01 
.  1  *613.646  ) 

• 19038358E+01 
.  19246693E+01 
.  19457222E  +  01 
. 1 9669990E+0 1 
.  1  9885045E  +  0 1 

•  201 024  39E  +  0  1 
. 2032222 1 t+0 1 
. 2 054444 6E+0 1 
•20769167E+01 

.  1 =649.684 ) 

.20996442E+01 
•21226330E+01 
.21 458891 E  +  0  1 
.21 694 1 88E+0 1 
.21 932286E+0 1 
• 22 1 73252 E+0 1 
.224 1 7 158E+0 1 

•  2266407SE+0 1 
. 229 1 4078E+0 1 

.  I  *685.  "*20  ) 

• 23167.. 46E+01 
• 23423660E+0 1 
•23683404E+0 1 

•  23946565E+0 1 
.242 1 3236E+0 1 

•  244835 1 0E+0 1 
•24757487E+01 
. 25035269E+0 1 
• 253 1 6964 E+0 1 

.  I *721.756) 

•  25602683E+0 1 
. 2589254 3E+01 
•261 B6666E+0 1 
.264851 79E+01 
•267882 16E+01 

•  2709591 4E+0 1 

•  27408420E  +  0 1 

•  27725887E  +  0  1 

•  280484 75E  +  0 1 
1 *757.792 ) 


l ( .28293877E+01 . 

2  .28625835E+0 1 . 

3  .28963395E+01 . 

4  . 2930673 1 E+0 1 . 

5  .296561 05E+01 . 

6  .3001 1670E+01 . 

7  . 3037367 1 E+0 1 . 


3  . 32269009E+0 1 • 

4  .326951146+01. 

5  .331096376+01. 

6  .335329336+01. 


.283763516+01 
.287096926+01 
. 2904868 3E +01 
.293935 1 9E  +  0 1 
.297444066+01 
•30101 SS8E+0 I 
.304652046+01 

.308355856+01 


.315975826+01 

.319897526+01 

.323897656+01 

.327979426+01 

.332146246+01 

.336401726+01 


. 1867891 3E+01 
. 18883519E+01 

. 1 9C90239E  +0 1 
. 192991 18E+01 
. 1 951 0202E+0 1 
. 19723537E+01 
. 19939173E+01 
• 20 1 571 59E+0 1 
. 20377546E  +  0 1 
•20600390E+01 
•20825744E+01 

•  2 1 053667E+0 1 
.21 2842 1 7E  +0 1 
.2151 7456E+01 

•  2 1 75344 7E  +0 1 
.21 992256E+0 1 
.222339516+01 
•22478602E+0 1 
. 22726263E+0 I 
.229770706+01 

•  2323 1 042E+0 1 
• 2348828 0E +01 
•23748870E+0 1 
. 2401 2900E+0 1 
•24280463E+01 
•24551653E+01 
. 24826572E+0 1 
•25105322E+01 
.253880126+01 

. 25674755E+0 1 
•25965670E+01 
• 26260878E+0 1 
. 26560509E  +0 1 
• 26864697E+0 1 
•271 73584E+0 1 
•2748731 6E+0 1 
•27806048E+01 
. 28 129941 E+0 1 

.284591676+01 
.287939036+01 
•291 34337E+01 
.294806666+01 
.298330986+01 
•30191852E+01 
.305571596+01 

•  30929262E  +0 1 
.313084216+01 

. 316949066+0 1 
.320890076+01 
.324910316+01 
•32901302E+0! 
.333201656+01 
.337479896+01 


. 1 8729869E+0 1 . 
. ie934999E+01 ) 

•  1 9 1 42255E  +  0 1 « 

•  1 935 1 68 1 E  +  0  1  . 
. 19563323E+01 . 
.197772286+01. 
. 1 9993447E+0 1 . 
.202120286+01 . 

•  20433025E  +  0 1  * 
• 2065649 1 E+0 1 . 

•  20882482E+0  1  ) 

•21 1 1 1056E+01 . 
•21 342272E+0 1 . 
.21 576 1 93E+0 1 . 
•2181 2882E+0 1 . 
•22052406E+01 . 
•22294833E+0 1 . 
.225402356+01. 
•22788686E+0 1 . 
.2304026 1 E+0 1 ) 

•23295042E+01. 
.235531 10E+01. 
.238 1 4552E+0 1 . 
•24079456E+0 1 • 
• 2434791 6E+0 1 . 
•24620030E+0 1 » 
•24895896E+0 1 . 
•251 7562 1 E+0 1 . 
•25459314E+0 1 > 

•  25747088E+0 1  * 

•  26039064E+0  1 « 
•26335366E+0 1 . 
•266361246+01. 
.2694 1 473E+0 1 • 
.2725 1 557E+0 1 . 
• 27566524E+0 1 • 
.278865316+01. 
.2821 1741E  +  01  ) 

.285423276+01. 
.288784696+01. 
•292203586+01. 
•295681936+01. 
.299221856+01. 
•302625556+01. 
.306495376+01. 
•310233806+01. 
.314043446+01  1 

.317927066+01 . 
•321887586+01. 
•325928126+01. 
.330051986+01 . 
.334262666+01. 

•  338563906+0  I « 


8  • 30742345E +0 1 . 

9  .311179436  +  01*  • 3 1 2  I 2955E  +  0 1 
DATA  I TABLEC 1 II . I *793,828 1 

1 (.31 500730E+0 1 
2  .318909866+01 


BC001780 
BC001790 
BC001600 
BC001810 
BC001820 
BC001830 
BC001840 
BC001850 
BC001860 
BC001870 
BC001880 
BC001890 
BC001900 
BC001910 
BC001920 
BC001930 
BC001940 
BC001950 
BC001960 
BC001970 
BC001980 
BC001990 
BC002000 
BC002010 
BC002020 
BC002030 
BC002040 
BC002050 
BC002060 
BC002070 
BC002080 
BC002090 
BC002100 
BC0021 10 
BC002120 
BC002130 
BC002140 
BC002150 
BC002160 
BC002170 
BC002I80 
BC002190 
BC002200 
BC0022I0 
BC002220 
BC002230 
BC002240 
BC002250 
BC002260 
BC002270 
BC002280 
BC002290 
BC002300 
BC00231 0 
8C002320 
BC002330 
BC002340 
8C0023S0 
60302360 
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7  .33965383E+01  .  .  34074972E+0  1  , 

8  .34407389E+01 .  • 345 1 9435E+0 l , 

9  .348S9386E+01 •  . 34974000E+0 1 . 
DATA  (TABLEC < IS • I ■829*864) 

11.353218346+01.  • 354391 37E+0 1 . 

2  • 35795229E+0 1 «  .3591 5350E+0 1 . 

3  . 36280  I 02E+0 1 «  . 36403 1 79E+0 1 . 

4  .36777022E+01 .  . 36903205E+0 1 . 

5  .37286603E+01 .  . 374 1 6054E+0 1 . 

6  • 37809509E+0 1  «  • 37942400E+0 1 . 

7  . 38346454E  +  0 1 .  . 38482973E+0 1 . 

8  . 388982 13E+01 .  . 39038564E+0 1 , 

9  » 39465627E+0 l •  • 3961 0032E+0 l < 
DATA  (TABLEC ( I  1 . 1 ■865.900 1 

1 (.40049610E+0! .  . 40 1 983 1 0E  +  0 1 , 

2  .40651  159E+01  .  . 408044 1 7E+0 1 « 

3  .4127! 364E+01 .  . 4 1 429467E+0 1 , 

4  .4191 141BE+01 .  • 42074685E+0 1 i 

5  .42572636E+0 1  «  .  4274  1  4  1  3E+0  1  i 

6  . 4 3236463E+0 1 <  . 4343 1 1 37E+0 1 , 

7  • 43964502E+0 1 .  .441 45498E+0 1 « 

8  .44698529E+0 1 «  .44886324E+01 < 

9  .4S460526E+01 .  . 45655649E+0 1 < 
DATA  (TABLEC ( 1 ). I >901 .936) 

1 < .462S2709E+01 .  .46455756E+01 < 

2  • 47077568E+0 1 .  . 472892 1 0E+0 1 , 

3  .4793791 5E  +  01 .  . 4B 1 389 1 2E+0 1 < 

4  .40836943E+O 1 .  . 49068 1 60E+0 1 , 

5  • 49778293E+0 1 .  . 5002072 1 E+0 1 , 

6  • 50766 1 49E+0 1 .  . 5 1 020929E+0 1 « 

7  • 51 805343E+0 1 .  . 52O73804E+0 1 . 

8  • 5290 1 5G0E+O 1  .  .  53 1  8520  1  E  +  0  1  • 

9  • 5406 1 253E+0 1 ♦  . 54 3620 1 1 E+0 1 
DATA  (  T  ABlEC ( I  1  . I »937 . 972 ) 

1 (.5529241 1E+01 .  . 556 1 24 1 8E+0 1 < 

2  • 5660 4337E+0 1 •  . 56946243E+0 I < 

3  . 580084 4 2E  +01.  . 58375425E+0 1 . 

4  • 595 1 8593E  *0 1 .  . 599 1 4643E+C 1 < 

5  .611 52 1 54E  +  0 1 .  .61582278E+01 

6  .62931 1 03E+01 .  . 6340 1 7 1 3E+0 1 

7  •  6488 3073E +  0 1 .  . 65403382E  +  C  1 

8  .6704(1 1  44E+CM  .  .67627B95E  +  01 

9  .69476761E+n;  ,  .70131  l 5PE  +  P 1 

DATA  ( TAHLFO (  I  1 .  I  ^9^3.999 ) 

1  (  .722  7  8'<6«>F.  ♦  *'  1  ,  .729931  76E.0  1 

2  .75445221L.C1  .  .  T63  il42‘i7E*C  1 

3  .79266326!  .*31  .  .  0  757  JO  0  1 

4  .83994  1  02F.  +  C  1  .  •<?*.  3739'  9E  +  01 

5  .  90  1  972C  C;-.  *0  1  .  .'-2I234C4E  +  0I 

(  .992  36  >  ' '•CO  1  .  .  1  T’|  99 je ♦  02 

7  .  1  1618206E*'.  .134^9716*  *  ?. 

EQuIvAl. NC!  (OPS.ICilf.l 

call  '-.ft  : : 

CALL  CRT  ID  (AHA39C) 

C  006  =  5 . /6. 

C6JF"6./F. 

L  A  =  1 
L0=  1 

5'  READ  1  _><.  .  J0B.FI4T 


. 34185165E+01 
.346321 1  I E+0 1 
. 35089274E+0 1 

.355571  31  E+01 
. 36036 196E+01 
. 3652701 8E+0 1 
. 37030 189E+01 
.37546347E+01 
•38076179E+01 
•  38620431 E  +  01 
.391 79908E+0I 
. 39755487E+0 1 

.40348123E+01 
.4095B857E+01 
.41 588831E+01 
•42239295E+0 1 
.4291  1627E+01 
•43607349E+01 
•44328 1 48E+0 1 
.45075899E+01 
• 45052695E+0 1 

•46660886E+0 1 
.475031 1 6E+0 1 
. 483B237BE+01 
.4 930208 0E+0 1 
.502661 22E  +  0 1 
• 5 1 278997E+0 1 
•5234591 7E+0 1 
.53472975E+01 
• 54667 360E+0 1 

• S5937628E+0 1 
• 5729406 0E+0 1 
•5B749267E+0 1 
. 6031 8700E+0  1 
.6202 1 856E+0 1 
.  6  388  366  4  E  +  0  1 
•6593 674 7E+0  1 
•68224954E+01 
.708091 B9E .01 

. 73777589E+01 
.  7726 4 6^76  +  01 
.81 490839E  +  0 1 
. 060  E  6  I  18E+01 
.942 1 061 4E  +0  I 
. 1 0596635F+02 
.1381551  IE +02 


. 3429596 9E+0  l.t 
.34745426E+01 . 
.3520521 6E+01  ) 

•3567S826E+01 . 
•36157777E+0  1 . 

•  36651629E  +  0  1  . 
. 37157985E+0 1 . 

•  37677495E+0  1  . 
. 382 1 0860E+0 1 . 
. 3875884 0E+01  . 
•39322257E+0 1 . 
.39902008E+0 1  ) 

.40499067E+01. 
.411 1 4500E+0 1 . 
.41 749474E+0 1 . 
.42405271 E+0 1 . 
•43083302E+0 1 . 
.43785 1 28E+0 1 . 
•44512481 E+0 1 . 
•4526728BE+0 1 « 
•46051 702E+0 1 ) 

.46868 142E+01 . 
• 477 1 9334E+0 1 . 
•4B608369E+0 1 . 
•4953B770E+0 1 . 
•50514573E+0 1 . 
•51540439E+0 I . 
.5262 1 733E+0 1 . 
.53764951E+01. 
. 54977444E+0 1 ) 

.5626821 4E+0 1 . 
•57648072E+0 1 • 

•  59 13023 1E+0 1  % 
.60731 0B5E+01. 

•  6247131 3E+0 1  * 
.6437751 6E+01. 
•664BA727E+0 1 . 

•  6884 0388E+0 1 . 
.9|51 1015E+01 > 

. 74594029E  +  0 1  * 
.78240460E+01 . 
•8270333 1 E+0 1 . 
•88456973E+0 1 . 
.O6566275E+0 1 . 
.1  1 04  2922E  +  02  » 


BC002370 
8C002380 
BC002390 
BC002400 
BCO02410 
BC002420 
BC002430 
BC002440 
BC002A50 
BC002460 
BC002470 
BC002480 
BC002490 
BC002500 
BC0025 1 0 

BC002520 
BC002530 
BC002540 
BC002550 
BC002560 
BC002570 
BC002580 
BC002590 
BC002600 
BC0026 1 0 
BC002620 
BC002630 
BC002640 
BC002650 
BC002660 
BC002670 
BC0026B0 
BC002690 
BC002700 
BC002710 
BC002720 
8C002730 
BC002740 
BC002750 
BC002760 
BC002770 
BC002780 
BC002790 
BC002800 
BC0028;0 
BC002820 
BC002B30 
BC0028v'0 
BC002850 
BC002860 
BC002870 
BC002S80 
BC002890 
BC  002900 
0COO291 0 
BC002920 
BC002930 

BC00294  0 

BC002950 
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BEST  AVAILABLE  COPY 


h  OQM  A  T  (  !  CAfl/K'AB  ) 
READ  1 1 5. NRUN 
FORMAT  (15) 

DO  SO  1  R'JN  -  1  « NRUN 

PRINT  110,  JOB  .FMT 


BC002960 

BC002970 

BC002980 

BC002990 

BC003000 


(  1H1 ,56X.  13HC0NTR0L  CARDS// 1 0A8//2X , 20HVAR I  ABLE  FORMAT  CARDBC003010 


1  ,5x  .  1  . .AS  ) 

PRINT  11 7, NRUN 

117  FORMAT  (IHO.  6HNRUN  =.13) 

120  RE  AC  15C.NPRINT. JPLOT.NOTIC.  INPUT  .NTYPES, ITRAN1  , ITRAN2, IRA, IRO* 
l  NP A  1 R .NN.NL.  (PS  <  JJ 1  , PDELT ( J J  > , PE  <  J J ) «  J J  = 1  , NN ) 

1 5"  FORMAT  < 1 2 1 5/1SF5.3 ) 

PRINT  155.1 RUN.NPR 1  NT , JPLOT .NOT  I C*  I NPUT , NT YPES , I TRAN1  ,  I TRAN2,  IRA 


I  RO  •  NP  A  I  R  ,NN,Nl_ 

(1H0.6HIRUN  = , I 3/1H0.9HNPR I  NT  =  . I l . 5X • 8H JPLOT  ■  ,IJ,5X. 

BHN0T1C  =  ,12, 5X , AH  1 NPUT  =  , I  1  , 5X , BHNT YPES  ■  , I1,5X, 

OH  I  T  RAN  1  =  ,  !2/lHO,9HlTRAN2  =  ,  I  2 , 4X  ,  6H I  flA  «  ,11, 7X. 

6HIR0  =  .11.  ex.SHNPAIR  =  .I2.4X.5HNN  «  ,11,9X.5HNL  ■  . I 

1 57, (PS ( JJ ) .PDELT ( JJ) .PE ( JJ ) . JJ»1 ,NN> 


BC003020 
BC003030 
BC003040 
3.  BC003050 

BC003060 
BC003070 
IRA.  BC0030BO 
BC003090 

X.  BC0031 00 
BC0031 10 
BC003I20 

.11)  BC0031 30 
BC0031 4Q 


157  FORMAT  (  1  HO . 43HP— VALUES  IDENTIFICATION  FOR  CHI-SQUARE  TEST/ 1H  ,5X,BC003150 
1  5(F5.3. IH < ,F5.3. 1H) .F5.3.2X ) )  BC003160 

IF  (JPLOT, LE.O)  GO  TO  163  BC003170 

IF  (JPLOT. FQ. 2)  GO  TO  163  BC003I80 

READ  160.  ( PPS ( J J ) , PPDELT ( J J ) , PPE ( J J ) , J J« 1 , NL  >  BC003190 

160  FORMAT  ((15F5.3))  BC003200 

PRINT  1 62, (PPS (JJ ) .PPDELT (JJ > .PPE (JJ ) ,JJ*I ,NL )  BC0Q3210 

162  FORMAT  ( 1  HO . 33HP— VALUES  IDENTIFICATION  FOR  PLOTS/IH  ,5X.5(F5.3.  BC003220 

1  |H(,F5.3«1H),F5.3,2X  ))  BC003230 

163  IF  (NTYPES.LE.O)  GO  TO  165  BC003240 

READ  164.  X1D1.XID2  BC003250 

164  FORMAT  (10A8/10A8)  BC003260 

•65  IF  ( ITRAN1 .EO.l )  GO  TO  168  BC003270 

READ  166.ATRAN1 .BTRAN1 .CTRANI .0TRAN1 .ETRAN1  CC003280 

166  FORMAT  (5E14.8)  BC003290 

PRINT  167, 1TRAN1 , A TRAN  1  , BTRAN1  .CTRANI .DTRANI »  ETRANI  BC003300 

167  FORMAT  ( 1  HO , 37HTRANSF0RMAT I  ON  CODE  FOR  XI  VALUES  IS  .12/IH  •  BC003310 

1  4  3HTRANSF  ORMAT I  ON  CONSTANTS  FOR  XI  VALUES  ARE— /IH  ,3X»  BC003320 

2  8HATRAN!  =,E15.8,3X. 8HBTRAN1  • . E 1 5. 8 . 3X . 8HCTRAN1  «,E15.8  BC003330 

3  , 3X , 8HDTRAN 1  a . E 1 5. 8 . 3X . BHETRAN1  ■•E15.H)  BC003340 

168  IF ( I TRAN2.EQ. 1 )  GO  TO  180  BC003350 

READ  1 69 , A  TRAN2 , 8TRAN2 , C  TRAN2 .DTRAN2 , ETRAN2  BCC03360 

169  FORMAT  (5E14.8)  BC003370 

PRINT  173,1 TRAN2 , ATR AN2 . BTRAN2 , CTRAN2 , DTRAN2 .ETPAN2  BC003380 

170  FORMAT  ( 1  HO , 37HTRANSF0RMAT I  ON  CODE  FOR  X2  VALUES  IS  »12/1H  ,  BC003390 

1  43H TRANSF ORMA T I  ON  CONSTANTS  FOR  X2  VALUES  ARE—/ 1 H  ,3X.  BC0034Q0 

2  8HATRAN2  = . E 1 5 . 8 . 3X . 8HB TRA N2  * , E l 5 . 8 , 3X .8HCTRAN2  >.E13.8«  BC003410 

3  3X.8HDTRAN2  = , E 1 5 • 8 . 3X , 8HETRAN2  -.E15.8)  BC003420 

1 80  PRINT  190  BC003430 

190  FORMAT  (  1H1 , 15X, 13H0RIGINAL  OAT A , 1 9X , 1 6HTR ANSFORMED  DAT  A/ 1 H  ,3X.  8C003440 

l  l(| .8X.2HX1 . 14X.2HX2 . S5X.2HX1 , 14X.2HX2)  BC003450 

IF ( IRUN.GT. 1 )  GO  TO  701  BC003460 

200  Il«I  BC003470 

I  2  =  NP A  I R  BC003480 

IF  ( INPUT. LE.O)  GO  TO  500  BC003490 

300  READ  (5.FMT)  I  END .  ( X 1  ( I  ) . X2 <  I  )  • I » I l . 1 2 >  BC003500 

IF  (1EMD.NE.0)  GO  TO  400  BC0035I0 

11»11+NPAIR  BC003520 

12=1 2+  NPA I R  BC003530 

GO  TO  3C0  BC003540 
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400 

500 


600 

700 

701 

705 

710 

715 

720 

725 

730 

735 

740 

745 

750 

755 

760 

765 
77  0 

775 

700 

785 

790 


795 

800 

805 

810 

815 

820 

82b 

830 

835 

84  0 
8«5 


85o 

855 

860 


NPTS=1 1+IEND-l 
GO  TO  700 

READ  FMT,  I  END. (XI  < I  )  ,X2< I ) » 1*1 t . 12) 

IF  (1END.NE.0)  GO  TO  600 

1  1  Ml  +  NP  A  I  R 

12=1 2+NPA I R 

GO  TO  300 

NPTS  =11+ 1  END- 1 

CALL  SEND  (NPTS.X1 .X10R1G) 

CALL  SEND  (NPTS.X2.X20RIG) 

IF  (ITRAN1.LE.1)  GO  TO  825 
GO  TO  <  825 »  705 .715. 725. 735. 74b. 75 5, 765 
DO  710  X«1 .NPTS 

XI  (K)  =  AL0G(X1  <K  )  ) 

GO  TO  825 

DO  720  X=  1  , NPTS 

XI  <K)*AL0G(AL0G(X1 (X))  > 

GO  TO  825 

DO  730  Xs  1 .NPTS 

XI (K ) = AlOG ( ATRAN1 +X 1 (<) ) 

GO  TO  825 
DO  740  K  = 1 . NPTS 

XI  (X )  =  ALOG( 8 TRAN 1 +ALOG (CTRANJ  +X 1 (Kill 
GO  TO  825 
DO  750  X  =  1  , NPTS 
XI  (X ) =SQRT (XI  (X  )  ) 

GO  TO  825 
DO  760  x  = 1  .NPTS 
XI  ( X  )  =  1 .O/Xl  (X) 

GO  TO  82b 

DO  770  K= 1 « NPT S 

XI (X)=l .0/(DTRANl+Xl (K) ) 

GO  TO  825 
DO  780  X  =  1  ,  NPTS 

XI  ( X ) =  A  T  AN (XI  (K)/SQRT( 1 .O-Xl  (X>«#2)> 
GO  TO  825 
DO  79u  X  =  1 .NPTS 

XI (K)=2.G*ATAN(SQRT(X1 (K J l/SORT ( 1 . 0-X 1 

GO  TO  825 

DO  800  X  = 1 .NPTS 

XI  (K)  sXl  (X  1/ETRAN1 

GO  TO  825 

DO  810  K  = 1 .NPTS 

XI (X ! =S 1 N (X 1 (X) > 

GO  TO  825 


BC003550 

BC003560 

BC003570 

BC003580 

BC00359C 

BC003600 

BC003610 

BC003620 

BC003630 

BC00364Q 

BC003650 

775. 785, 795.805, 815) « ITRAN1BC003660 

BC003670 
BC003660 
BC0C3690 
BC003700 
BC00371 0 
BC003720 
BC003730 
BC003740 
BC003750 
BC003760 
BC003770 
BC003780 
8C003790 
8C003800 
BC00381 0 
BC003820 
BC003830 
BC003S40 
BC003850 
BC003860 
BC003870 
BC003880 
BC003890 
BC003900 
BC00391 0 

<>>>  BC003920 

BC003930 

BC003940 

BC003950 

3C003960 

BC003970 

BC003980 

BC003990 


DO  82o  X  =  l  .NPTS 
XI  ( X )  =  COS ( X 1  (X) ) 

IF  ( I TRAN2 »LE • 1  )  GO  TO  950 

GO  TO  (950, 830. 940, 850. 860. 870. 880,890.900. 91 
DO  835  K= 1 .NPTS 
X2(X)=AL0G(X2(K) ) 

GO  TO  950 

DO  845  K= 1 , NPTS 

X2 ( X  )  =  ALOG ( AlOG ( X2 (X  >  )  ) 

GO  TO  95C 

ro  855  x  =  1  . NPTS 

X2 (x )  =  AL0o ( ATRAN2  +  X2 (X  > ) 

GO  TO  950 
DO  865  X  =  1  , NP T S 


BC004000 
BC00401 0 
BC004020 

0,920,930.940 ) » I TRAN2BC004030 

BC004040 
BC004050 
BC004060 
BC004070 
8C004080 
BC00A090 
0COOA1OO 
BC004! 10 
BC004120 
BC004130 


53 


B65 

X2(K  )  *A|_OG (STRAN2+AL0G (CTRAN2+X2 (K  ) ) 1 

BC004140 

GO  TO  95G 

3C004150 

870 

DO  875  K  =  1 « NPTS 

BC004160 

875 

X2 ( K ) =  SORT ( X2 (K  )  ) 

BC004170 

GO  TO  950 

BC004180 

eeo 

DO  885  K=  1  <  NIPT 3 

BC004190 

885 

X2(K)=1 ,0/X2(K) 

BC004200 

GO  TO  950 

BC004210 

890 

DO  895  K  =  1 « NPTS 

BC004220 

895 

X2<K}=1 • 0/ ( DTRAN2+X2 (K ) ) 

BC004230 

GO  TO  950 

BC004240 

900 

DO  905  K=1  < NPTS 

BC004250 

905 

X2 ( K 1  * AT AN ( X2 <K  J/SQRT (1.0- 

-X2 (K 1**2  >  > 

BC004260 

GO  TO  950 

BC004270 

9!C 

DO  915  K-l » NPTS 

BC004280 

915 

X2<K)=2«G*AT'N( SORT ( X2 (K ) )/SQRT( 1 ,0-X2(K> ) ) 

BC004290 

GO  TO  950 

BC004300 

920 

DO  925  K=1,NPTS 

BC004310 

925 

X2 (K )=X2 (K 1/ETRAN2 

PC004320 

GO  TO  950 

BC004330 

930 

DO  935  K=1.NPTS 

BC004340 

935 

X2 (K ) *S IN ( X2 (K  )  1 

BC004350 

GO  TO  950 

SC004360 

940 

DO  945  K=1 « NPTS 

BCQ04370 

945 

X2(K)*C0S(X2(K) ) 

BC004380 

950 

1  PAGE  *50 

BC004390 

DO  960  1*1, NPTS 

BC004400 

IF  ( I  PAGE • NE .0 1  GO  TO  954 

BC0044 1 0 

I PAGE=50 

BC004420 

PRINT  951 

BC004430 

951 

FORMAT  (1H1) 

BC004440 

954 

PRINT  955,1  , X10R !G( I ) ,X20RIG( I ) , XI  < I  1  , X2 ( I  ) 

BC004450 

955 

FORMAT  (1H  , I4,2X,E14.8,2X,E14,8,3XiE14.8«2X,E14.B) 

BC004460 

960 

I  PAGE* I  PAGE- 1 

BC004470 

X1MIN=X1 ( 1 ) 

BC004480 

X2MIN=X2( 1 ) 

BC004490 

S I GX 1 2=X 1 < 1 1**2 

BC004500 

S I GX22=X2 < 1 1**2 

BC004510 

R*X 1 < 1  5*X2< l  ) 

BC00452P 

XIMEAN*X1 ( 1 ) 

BC004530 

X2MEAN*X2 ( 1 ) 

BC004540 

XI MAX*X 1 M I N 

BC004550 

X2MAX*X2M I N 

BC 004560 

DO  1350  I *2, NPTS 

BC004570 

IF  (XI ( I 1.GT.X1MAX)  GO  TO 

970 

BC004580 

IF  (XI <1 l.LE.XIMINl  GO  TO 

990 

BC004590 

GO  TO  1000 

BC004600 

970 

X1MAX-X1 ( I ) 

BC004610 

GO  TO  1000 

BC004620 

990 

XI M I N*X 1 ( I > 

BC 004630 

1000 

IF  (X2( 1 1.GT.X2MAX)  GO  TO 

1100 

BC00464Q 

IF  (X2 ( I1.LE.X2MIN)  GO  TO 

1200 

BC004650 

GO  TO  1300 

BC004660 

l  100 

X2MAX*X2( I 1 

BC004670 

GO  TO  1300 

BC0046S0 

1200 

X2M1N*X2 ( I 1 

BC004690 

1300 

S!GX12*SIGX12+X1 ( l )**2 

BC004700 

SIGX22«S I GX22+X2 ( 2 1**2 

BC004710 

R-R+Xl ( 1  )*X2 (  1  ) 

BC004720 
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X 1 MEAN*x 1 MEAN+X 1  (  1  )  BC004730 

X2M£AN»X2MEAN+X2 ( I >  BC004740 

1350  CONTINUE  BCQ04750 

PRINT  1250*  X1MIN.X1MAX.X2MIN.X2MAX  BCQ04760 

1250  FORMAT  ( 1 H 1 «8HMI N  XI  =  .£ 1 4 . 8 . 1  OX. 8HMAX  XI  «.E14.8/1H  .8HMIN  X2  =  ,  BC004770 
1  E14.8.10X.8HMAX  X2  ».E14,8)  BC004780 

XNPTS=NPTS  BC004790 

S I GX 1 2=S I GX 1 2-X 1 MEAN*»2/XNPTS  BC004800 

S I GX22=S 1 GX22-X2ME AN*#2/XNPTS  BC004810 

Rs (R*XNPTS-X1M£AN*X2MEAN)/(XNPTS*SQRTCSIGX12*SIGX22>  >  BC004820 

S I GX 1 2*  SIGX12/IXNPTS-1 . >  BC004830 

SI GX22=  S I GX22/ I XNPTS- 1 < )  BC004840 

X 1  MEAN  =  X 1 MEAN/XNPTS  BC004850 

X2MEAN  «  X2 MEAN/XNPTS  BC004860 

SIGX1  =  SORT ( S I GX 1 2 )  BC004870 

SIGX2  «  SORT I S I GX22 )  BC004880 

PRINT  1400.R.XI MEAN . X2ME  AN. S I GX 1 2  « S I GX22  *  S I GX 1 .SI  6X2  BCo0489o 

1400  FORMAT  ( 1H0.2HRs.E14. 8. 3X.8HMEAN  X 1 = . E 1 4 .8 . 4X . 8HMEAN  X2» . E 14. 8, 4X. BC004900 

1  1 2HVARI ANCE  X 1 = . E 1 4 .8 ,4X ♦ 1 2HVAR I ANCE  X2».E14,8./  1H  ,  BC004910 

2  12HSTAN  DEV  XI =.E14«6*5X. 12HSTAN  DEV  X2=.El4.8>  BC004920 

A1*R#SIGX2/SIGX1  BC004930 

A2=X2MEAN-A 1 #X1 MEAN  BC004940 

A3=R*S I GX 1 /S I GX2  BC004950 

A4  =  X 1  ME AN- A3 #X2 ME AN  BC004960 

PRINT  1450. A2.A1 .A4. A3  BC004970 

1450  FORMAT  I 1  HO , 3HX2= . E 1 4.8 « 3H  +  .E14.8.4H  *X1/4H  X1*.E14.8.3H  +  , E 1 4 .BC004980 
18. 4H  *X2 )  BC004990 

PRINT  1475  BC005000 

1475  FORMAT  <1H1)  BC005010 

NT  =  0  BC005020 

DO  1600  KN=1,NN  BC005030 

NT  =  NT  + 1  BC005040 

IF  (PDELT(KN).EO.O. )  GO  TO  1550  BC005050 

NT  1  *  (PE(KN)-PS(KN) )/PDELT(KN)+0.0001  BC005060 


PINT ) =  PSIKN) 

NTP 1  rNT+1 
NT  «  NT+NT1 
PRO0  «  PS(KN) 

DO  1500  KM*NTP1 .NT 
PROB  *  PROB+PDELT (KN) 

1500  P ( KM ) *  PROB 
GO  TO  1600 
1550  P I  NT ) =  PS (KN ) 

1600  CONTINUE 
NPROB  =NT 
DO  1700  LMa 1 . NT 
IS  ■  P (LM >*1 000.+.0005 
1700  C2(cM+l)  ■  TABLECIIS) 

NT«NT  + 1 

RCOMPs 1 ./ ( 1 .-R**2 ) 

DO  1800  Isl.NPTS 

AA=(  (Xt  (  I  ) -X 1  MEAN ) /S I GX 1  ) 

AB«(  1X2  I  I  )-X2MEAN)/SIGX2 ) 

1800  0(1)  »RCOMP*  (  4A*#2-2,*R-»AA*AB  +  A8*»2  > 
C2 I  1  )«C 

IF  INPRINT. LE.O)  GO  TO  2650 
I L t  NE  »50 
DO  2200  J»2,NT 
JMI ■ J-l 


BC005070 
BC005080 
BC005090 
BC005100 
BC0051 10 
BC005120 
BC005130 
BC005140 
BC005 1 50 
BC005160 
BC005I 70 
BC005I80 
BC005190 
BC005200 
BC005210 
BC005220 
BC005230 
BC0052«  0 
BC005250 
BC005260 
BC005270 
BC005280 
BC005290 
BC005300 
BC005310 


1850 
1  875 
1  900 

20C0 


2050 


21  CO 


22C0 


2250 

2300 


2400 


2500 

2650 

2700 


2750 


2800 


2850 


2900 

2975 


3000 

3100 


15  (ILINE.NE.O)  GO  TO  1875  6C005320 

I L I NE=50  BC005330 

PRINT  1850  BC005340 

FORMAT  ( 1 H 1 )  BC005350 

PRINT  1900, JM1 ,C2(J).JM1 ,P( JM1 )  BC005360 

FORMAT  (1H  ,4H  C2(,13,4H)  *  , E 1 4. 8 . 5X 1 2HP ( . I  3 , 4H )  *  ,F8.3)  BC005370 

I L I NE  = I L I NE- l  BC005380 

DO  2200  1=1* NPTS  BC005390 

IF  (ILINE.NE.O)  GO  TO  2050  BC005400 

I L I NE=50  BC005410 

PRINT  1850  BC005420 

IF  ( (C2(J-1 ).GT.Q< 1 ) ).0R  .  < C2 < J ) .LT .Q ( 1 ) > >  GO  TO  2200  BC005430 

PRINT  2100,1  *X1  ( I  ) • X2 < I  )  »Q(  I  )  BC005440 

FORMAT  (1H  , !5,2H-(,E14.8.1H, . E 1 4 .8 » 2H ) -, E 1 4 .8 )  BC005450 

IL1NE= ILINE-1  BC005460 

CONTINUE  BC005470 

IF  (ILINE.NE.O)  GO  TO  2250  BC005480 

I L I NE=50  BC005490 

PRINT  1850  BC005500 

PRINT  2300. NT, NT  BC005510 

FORMAT  (IH0.4M  C2(,I3,12H)  *  INF  I N I T Y . 1 1 X . 2HP< . 1 3, 1 2H )  *  1 . 000 1BC005520 

IL I NE=1 LINE-1  BC005530 

DO  2500  1*1, NPTS  BC005540 

IF  (ILINE.NE.O)  GO  TO  2400  BC005550 

IL!NE=50  BC005560 

PRINT  1850  BC005570 

IF  (0(1) .LE.C2 (NT ) )  GO  TO  2500  BC005580 

PRINT  2100, I » XI  (I )  ,X2d  >  »Q( I )  BC005590 

I  LINE* I  LINE-1  BC005600 

CONTINUE  BC005610 

PRINT  2700. IRUN. NPTS, NT  BC005620 

FORMAT  (IH1 ,39x. 1 5HCH I -SQUARE  TEST/1  HO. 18X« 1 1HRUN  NUMBER  .13.  BC005630 

1  6H  WITH  ,15,1 8H  OBSERVATIONS  AND  .I4.10H-  INTERVALS/1H0,  BC005640 

2  8HINTERVAL.3X, 1 1HUPPER  B0UND.6X ♦ 1 1 HPROBAB IL 1  TV ,3X.  BC005650 

3  8  HOBS  FREO , 7X , 9HTHE0  FREQ.6X. 13HCHJ-SQU  CONTR  /1H  BC005660 

4  8  H I NDEX ( J  > , 6X , 5HC2 ( J ) , 13X  »4HP ( J )/ )  BC005670 

P (NT) =1.000  BC005680 

CALL  MOVE  (0, IOBS( 1 ) ,999)  BC005690 

00  2850  I  P  =  1  . NPTS  BC005700 

DO  2800  I C  *2 ,NT  BC005710 

IF  (C2( ICj.GT.Ot |P),AND.C2( IC-1 l.LE.O(IP) 100  TO  2750  BC005720 

GO  TO  2800  BC005730 

IOBS(IC-) )-IOBS(IC-l )+l  BC005740 

GO  TO  2850  BC005750 

CONTINUE  BC005760 

IOBS(Nt)  •  I0BS<NT)+1  BC005770 

CONTINUE  BC0057B0 

E(1 )«XNPTS»P(I i  BC005790 

DO  2900  IE *2, NT  BC005B00 

EIIEI  *  XNPTS*  (PdE»-P(  IE-1  »  )  BC005S10 

CSUM*0  BC005820 

DO  3000  I J* I »NT  BC005830 

IC  *IJ+1  BC00S840 

CSC  «<OBSd J)-E(IJ) )«*2/E(IJ)  BC005850 

CSUM  ■  CSUM  +CSC  BC 009880 

PRINT  3100.  IJ.C2(IC).P( I J) . IOBS( IJI.EIIJI .CSC  BC00SB70 

FORMAT  (1H  . 2X, I  3 ,5X «E 14 *8«SX «F8. 3, 7X« 1 4 . 1 0X.F6.4 .9X.F6.4 )  BC009B80 

K I -NT-6  BC009B90 

IF(KI.LE.O)  GO  TO  3710  BC009900 
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PRIM  3700 * C SUM  »K  1 

3700  FORMAT  (  I  HO  ,  1 2HCH I -SQUARE  «.F12.6.6M  WITH 
I0M/1H2 ) 

IF(Kl.GT.O)  GO  TO  3730 
3710  PRINT  3725 


BC00S910 

I3.19H  DEGREES  OF  FREEDBC005920 

BC005930 

BC005940 

BC005950 


3725  FORMAT  ( 1  HO • 76HCH 1 -SOUARE  COULD  NOT  BE  COMPUTED  BECAUSE  OF  INSUFF IBC005960 
1 C I  ENT  NUMBER  OF  INTERVALS)  BC005970 

3730  IF  (JPLOT.LE.O)  GO  TO  9000  BC0059BO 

IF  (JPLOT.EQ.2)  GO  TO  3826  BC005990 

NP*0  BC006000 

DO  3800  KN*1.NL  BC006010 

NP  =  NP  + 1  BC006020 

1 F  (PPDELT (KN) .EQ.O. )  GO  TO  3775  BC006030 

NP  I  ■  ( PPE  (KN )  -PPS  (KN )  )  /PPDELT  I KN  >4-0*0001  BC006040 

PP(NP)=PPS(KN)  BC006050 

NPPi  =  NP+ 1  BC006060 

NP  *  NP+NP1  BC006070 

PROB1  *  PPS (KN )  BC006080 

DO  3750  «M«NPP1.NP  BC0C6090 

PROS  I  *  PR081 +PPDELT  (KN )  BC006100 

3750  PP(KM)  =PROB 1  BC006110 

GO  TO  3800  BC006120 

3775  PP ( NP ) =  PPS (KN )  BC006130 

3800  CONTINUE  BC006140 

NPROBI  *  NP  BC006150 

DO  3825  LP*1.NP  BC006160 

ISP*PP(LP)»1000.+«0005  BC006170 

3825  C2P ( LP ) “TABLEC ( I SP )  BC006180 

GO  TO  3827  BC006I90 

3826  IF  (NTYPES.EQ.I >  GO  TO  3829  BC006200 

PRINT  3828. JOB, IRUN.NPTS  8C006210 

3828  FORMAT (2HS2, 10A8/2HS  .1IHRUN  NUMBER  .I3.6H  WITH  .I4.13H  OBSERVAT I OBC006220 

INS)  BC006230 

GO  TO  4100  BC006240 

3827  IF  (NTYPES.LE.O)  GO  TO  3840  BC006250 

3829  PRINT  3830.XID1 .XID2  BC006260 

3830  FORMAT  I 2H*2 . 1 0A8/2HS  « 1 0A8 )  BC006270 

GO  TO  4100  BC006280 

3840  PRINT  3850. JOB, 1  RUN, NPTS . NP. (PPS < JJ > .PPDELT ( JJ >, PPE ( JJ >» JJ- 1 . NL  BC006290 
1  )  BC006300 

3850  FORMAT  < 2H*2 , I 0AB/2H*  ,11 HRUN  NUMBER  .I3.6H  WITH  , I  4 , 1 8H  OBSERVAT I BCQ063 1 0 


IONS  AND  .I3.9H  CON TOUR S/2HS  « 5 ( F5 .3. 1 H ( ,F5. 3. I H > ,F5. 3 ) ) 
4J00  IF  (NOTIC.LE.O)  GO  TO  4200 
PNOTIC  *  NOTIC 
GO  TO  4300 
4200  PNOTIC  «  15. 

4 3C0  XIRANG  =X1MAX-X|MIN 
X2RANG  *X2MAX-X2MIN 
R 1 e  AMA  x 1  ( X 1  PANG , X2RANG ) 

D«  R1 /PNOTIC 

IF  (O.LT.l . )  GO  TO  4900 

DO  4500  IP-1 .9 

IF  (O.EQ.TABPt IP) >G0  TO  4600 

IF  (  T  Ag)P  (  I P  )  .LT.D.  AND.  TA8P  <1P+I)«GT.D)  GO  TO  4700 
4507  CONTINUE 

IF  (D.GT.TABPIIO) )  D-TABP(IO) 

GO  TO  5300 

4603  I D“T  A8P ( I P )  +  , 005 
GO  TO  4750 


BC006320 

BC006330 

BC006340 

BC006350 

BC006360 

BC006370 

BC006380 

BC006390 

BC006400 

BC006410 

8C006420 

BC006430 

BC006440 

BC006450 

BC006460 

BC006470 

BC006480 

BC006490 


4700  IDI «D/TABP< IP)+.9999 

D»TABP( I P ) *FL0AT ( 101 > 

GO  TO  5300 
4750  0“ I D 

GO  TO  5300 

4900  IF  (D.LT.TABNII ) )  GO  TO  5100 
DO  5000  I NEG« 1*9 
IF  (D.EQ.TABN! INEG) )G0  TO  5100 

IF  (TABNI INEGl.LT.O.AND.TABNI INEG+1 ) .GT.D )G0  TO  5200 
5000  CONTINUE 
5100  D« TABNI 1 ) 

GO  TO  5300 

5200  1D2«D/TABNI INEG )  +.9999 

D  «TABN( INEG 1*FL0AT< 1021 
5300  DX*D 
DY  *0 

X0R-X1MIN 

Y0R-X2MIN 

3875  YP*C5D6*  (XI MAX-X 1 M  l N )+X2MI N 
IF I YP— X2MAX )  4000.5375.3900 
3900  X2MAX*YP 

GO  TO  5375 

4000  X1MAX«C6D5#  ( X2MAX-X2M IN )+Xl MI N 

5375  CALL  GRF  (XI M IN .X2M IN.X 1  MAX .X2MAX.X0R. YOR. DX .OY . IRA. IRO.O.NPTS. 

1  XI (1 ).X2(1).LA,L0.2H54,-1.,0> 

IF  (JPL0T.EQ.2)  GO  TO  9000 
IF  (JPLOT.LE.O)  GO  TO  9000 
ORT 1*1X1  MAX— XI M IN }/4. 

XU  *X1MIN+QRT1 
XI 2  -X1MAX-ORT1 
HANOY (5 )*X1 MEAN 
HANDY  16 J-X2MEAN 

TEMPI  «  (X 1 1-X1 ME AN )*S I GX2/S I GX 1 

TEMP2  b ( X 1 2-X 1  ME AN ) *S I GX2/S I GX 1 

HANDY(2)*X2MEAN+R  *TEMP1 

HANDY ( 4 ) -X2MEAN+TEMP1 /R 

HANDY  IB )«X2MEAN+R  *TEMP2 

HANDY ( l 0 1 «X2MEAN+TEMP2/R 

HANDY! 1 )*X1 1 

HANDY (3)«X1 1 

HANDY  ( 7 ) «X 1 2 

HANDY (9) *X12 

0RT2« (X2MAX-X2MIN)/4, 

X21 "X2MIN+ORT2 
X22»X2MAX-QRT2 
HANDY! 12) aX2l 
HANDY! 14)*X21 
HANDY (16) *X22 
HANDY! 18 )«X22 

TEMP3» (X2I-X2MEAN)*SIGX1/SIGX2 
TEMPA« (X22-X2MEAN J#SIGXI/SIGX2 
HANDY! 1 1 >»XIMEAN+R  4TEMP3 
HANOY ( 1 3 ) -X l MEAN+TEMP3/R 
HANDY! 15) aXlMEAN+R  * TEMP 4 
HANOY! 17)»XIMEAN+TEMP4/R 

CALL  GRF  (0..0..0..0..0..0..0..0. .0.0.1 .9* HANOY ( I  ) .HANDY <2 ) .2*2* 
1  2H00.-1.0) 

X0MR2»I.-R  *#2 
DO  5500  Ja | ,WP 


BC006S00 

BC006510 

BC0C5520 

BC006530 

BC006540 

BC006550 

BC006560 

BC006570 

BC006580 

BC006590 

BC006600 

BC006QH0 

BC006620 

BC006630 

BC006640 

BC006650 

BC006660 

BC006670 

BC006680 

BC006690 

BC006700 

BC00671 0 

BC006720 

BC006730 

BC006740 

BC006750 

BC006760 

BC006770 

BC006780 

BC006790 

BC006800 

BC006810 

BC006820 

BC006830 

BC006840 

BC006850 

BC006860 

BC006870 

BC006880 

BC006890 

BC006900 

BC00691 0 

BC006920 

BC006930 

BC006940 

BC006950 

BC006960 

BC 006970 

BC006980 

BC006990 

BC007000 

BC00701 0 

BC 007020 

BC007030 

BC007040 

BC0070S0 

BC007060 

BC007070 

BC0070B0 
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V-C2PIJ)  BC007090 
TEMP 1  *  SIGX1 *SORT (V )  BC007100 
XINTVL-1 »999#TEMPI/499t  BC007I10 
X0R0-X1 WEAN-TEMPI -XI NT  VL  BC007120 


00  5400  L»l«500 

XORD*XORDTX I NTVL 

XTMP»  (  X0R0-X1 MEAN  >,  _,I6X1 

DISC* SORT ( ABS ( < V-XTMP442 )*X0MR2  J ) 

XI <2#L  >«XORD 

XI <2*L-1 > ■XORD 


BC007130 

BC007140 

BC007I50 

BC007160 

BC007I70 

BC0071B0 


X2<2#L-1 >»x2m£AN+SIGx2*(R  4XTMP+0 I SC ) 

5400  X2 (2#L ) “X2MEAN+S IGX2* (R  •XTMP-DISO 

CALL  GRF  (0».0..0..0..0.,0..0..0*«0*0«1  .1000**1 <1 >«X2<  »  1.1.1.2H75 

1  -1.0) 

5500  CONTINUE 
9000  CALL  I NTVL (V) 

PRINT  1 1 000* Y 

11000  FORMAT  (27H  RUNNING  TIME  IN  SECONDS  »  .E14.8) 

CALL  SEND<NPTS.X10RIG.X1  ) 


10000 


call  send 

CONTINUE 
GO  TO  *50 
STOP 
END 


(NPTS.X20RIG.X2) 


BC007200 
.BC007210 
BC007220 
BC007230 
BC007240 
BC007250 
BC007260 
BC007270 
BC007280 
BC 00 7290 
BC007300 
BC007310 
BC007320 
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'*■  IPONOOPINO  MILITARY  ACTIVITY 


A  description  is  Riven  of  the  EI-CIII  (Bivariate  Chi-Souare  Test)  computer 
program,  The  nrogram  is  for  assessing  the  agreement  oT*a  continuous  two- 
dimensional  sample  distribution  with  a  parent  bivariate  normal  distribution, 
fhe  program  computes  the  Chi-snuare  statistic  for  testing  the  null  hypothesis 
The  random  sample  is  from  a  bivariate  normal  parent  population."  The  minimum 
gnd  maximum  sample  values  of  each  random  variable  are  determined,  and  the  samrle 
estimates  of  the  five  parameters  (two  means,  two  variances,  and  the  correlation 
Coefficient)  of  the  bivariate  normal  distribution  are  computed.  Also,  the 
sample  estimates  of  the  two  regression  lines  (xj  on  xj  and  xj  on  xi)  arc  computed. 
An  optional  feature  allows  additional  output  of*a  plot  of  the  sample  data  and 
ghc  contour  ellipsos.  In  addition  to  an  example  problem  illustrating  the  input 
and  the  output  format,  a  flow  chart  and  program  listing  are  included  in  the  renort. 

The  program  is  coded  in  FORTRAN  IV  for  the  IBM  7030  (STRETCH)  committr. 
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